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Abstract 

The two-body hadronic decays of B mesons into pseudoscalar and axial vector mesons are studied 



C3 



within the framework of QCD factorization. The light-cone distribution amplitudes (LCDAs) for 



3 Pi and l P\ axial-vector mesons have been evaluated using the QCD sum rule method. Owing to 
the G-parity, the chiral-even two-parton light-cone distribution amplitudes of the 3 Pi Q~P\) mesons 
£T) ', are symmetric (antisymmetric) under the exchange of quark and anti-quark momentum fractions in 

the SU(3) limit. For chiral-odd LCDAs, it is other way around. The main results are the following: 
(i) The predicted rates for af (1260)^, 6^(1235)^^, 6?(1235)tt~, afK~ and bfK~ modes are in 
good agreement with the data. However, the naively expected ratios B(B~ — > a®n~) /B(B° — > 
afO < 1, B(B- -» a{7r°)/B(B° -» a^vr+) ~ \ and B{B~ -> b° 1 K-)/B{B° -» bfK~) ~ \ are 
not borne out by experiment. This should be clarified by the improved measurements of these 
decays, (ii) Since the B — > b\K decays receive sizable annihilation contributions, their rates 
are sensitive to the interference between penguin and annihilation terms. The measurement of 
B(B° — ► bfK~) implies a destructive interference which in turn indicates that the form factors 
for B — > b\ and B — » a\ transitions are of opposite signs, (iii) Sizable power corrections such as 
weak annihilation are needed to account for the observed rates of the penguin-dominated modes 
K{ (1270)vr+ and (1400)vr+. (iv) The decays B -> K{K with K x = K x (1270), K x (1400) are 
in general quite suppressed, of order 10~ 7 ~ 10~ 8 , except for B° — > Ki(1270)K° which can have 
a branching ratio of order 2.3 x 10~ 6 . The decay modes K^~K + and Kf K~ are of particular 
interest as they proceed only through weak annihilation, (v) The mixing-induced parameter S 
is predicted to be negative in the decays B° — > afir^, while it is positive experimentally. This 
may call for a larger unitarity angle 7 > 80°. (vi) Branching ratios for the decays B — > fiir, f±K, 
hnt and hiK with f x = /i(1285), /1 (1420) and h x = /ii(1170), h x (1380) are generally of order 
10~ 6 except for the color-suppressed modes f%iT and /ii7T° which are suppressed by one to two 
orders of magnitude. Measurements of the ratios B(B~ — > hi(1380)ir~)/B(B~ — > /ii(1170)-7r _ ) 
and B(B -> fx(U20)K)/B(B -> /i(1285)Z) will help determine the mixing angles 6i Pl and 9 3Pl , 
respectively. 
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I. INTRODUCTION 



In the quark model, two nonets of J p = 1 + axial-vector mesons are expected as the orbital 
excitation of the qq system. In terms of the spectroscopic notation 2S+1 Lj, there are two types 
of p-wave mesons, namely, 3 Pi and l P\. These two nonets have distinctive C quantum numbers, 
C = + and C = — , respectively. Experimentally, the J PC = 1 ++ nonet consists of ai(1260), 
/i(1285), /i(1420) and K 1A , while the l+~ nonet has &i(1235), /ii(1170), /ii(1380) and K 1B . The 
physical mass eigenstates K \ (1270) and ifi(1400) are a mixture of K\a and K\b states owing to 
the mass difference of the strange and non-strange light quarks. 

The production of the axial-vector mesons has been seen in the two-body hadronic D decays: 
D -» X"oi(1260), D° -»• iff (1270)vr and D+ -» iff (1400)tt , and in charmful B decays: B -» 
J/^iTi(1270) and B -» Doi(1260) [l]. As for charmless hadronic B decays, B° -> af (1260)tt t are 
the first modes measured by both B factories, BaBar and Belle. The BaBar result is 

B(B° -» af (1260)7r T ) = (33.2 ± 3.8 ± 3.0) x 10~ 6 . (1.1) 

where the assumption of B(af — ► 7r ± 7r ± 7r^) = 1/2 has been made. The Belle measurement gives 

B 

B(B° -» af (1260)tt t ) = (29.8 ± 3.2 ± 4.6) x 10~ 6 . (1.2) 

The average of the two experiments is 

B(B° -» af (1260)vr T ) = (31.7 ± 3.7) x 10~ 6 . (1.3) 

Moreover, BaBar has also measured the time-dependent CP asymmetries in B° — * a^(1260)7r T 
decays [4|. From the measured CP parameters, one can determine the decay rates of afir~ 
and a^7r + separately 0]. Recently, BaBar has reported the observation of the decays B° — » 
bfit^,bfK~ and B~ — > b\ir~ , b\K~ , o^vr - , aj~7r° 0, 0|. The preliminary BaBar results for 
B° -» K{ (1270>7r+ Kr(U00)ir+,afK-, B~ -> a7if°, /i(1285)if~, /i(1420)if- are also avail- 
able recently BflflQ- 

In the present work we will focus on the B decays involving an axial-vector meson A and a 
pseudoscalar meson P in the final state. Since the 3 Pi meson behaves similarly to the vector 
meson, it is naively expected that AP modes have similar rates as VP ones, for example, B(B° — ► 
a] t (1260)7r =F ) ~ B(B° — > p^ir^). However, this will not be the case for the l P\ meson. First of 
all, its decay constant vanishes in the SU(3) limit. For example, the decay constant vanishes for 
the neutral 65(1235) and is very small for the charged 6i(1235) states. This feature can be checked 
experimentally by measuring B° — > bfir~, 6^7r + decays and seeing if the former is suppressed 
relative to the latter. Second, its chiral-even two-parton light-cone distribution amplitude (LCDA) 
is anti-symmetric under the exchange of quark and anti-quark momentum fractions in the SU(3) 
limit due to the G parity, contrary to the symmetric behavior for the 3 Pi meson. 

Charmless B — > AP and B — > AV decays have been studied in the literature 
E3- Except for 0,0 , most of the existing calculations were carried out in the framework of either 
naive factorization or generalized factorization in which the nonfactorizable effects are described by 
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the parameter , the effective number of colors. In the approach of QCD factorization, nonfactor- 
izable effects such as vertex corrections, hard spectator interactions and annihilation contributions 
are calculable and have been considered in [ill . [l^ | for the decays B — > ai(1260)-7r, ai(1260)#, 
B -> fa(1235)iT\ &i(1235)#* and 6i(1235)p. 

One crucial ingredient in QCDF calculations is the LCDAs for 3 Pi and 1 P± axial- vector mesons. 
In general, the LCDAs are expressed in terms of the expansion of Gegenbauer moments which 
lave been systematically studied by one of us (K.C.Y.) using the light-cone sum rule method 
US]- Armed with the LCDAs, one is able to explore the nonfactorizable corrections to the naive 
factorization. 

The present paper is organized as follows. In Sec. II we summarize all the input parameters 
relevant to the present work, such as the mixing angles, decay constants, form factors and light-cone 
distribution amplitudes for 3 Pi and l P\ axial-vector mesons. We then apply QCD factorization 
in Sec. Ill to study B — » AP decays. Results and discussions are presented in Sec. IV. Sec. V 
contains our conclusions. The factorizable amplitudes of various B — * AV decays are summarized 
in Appendix A. 



II. INPUT PARAMETERS 
A. Mixing angles 

In the quark model, there are two different types of light axial vector mesons: 3 Pi and ^i, 
which carry the quantum numbers J PC = 1 ++ and l" 1 , respectively. The 1 ++ nonet consists of 
ai(1260), /i(1285), /i(1420) and K 1A , while the 1+" nonet has 61 (1235), &i(1170), fa (1380) and 
Kib- The non-strange axial vector mesons, for example, the neutral ai(1260) and 6i(1235) cannot 
have mixing because of the opposite C-parities. On the contrary, the strange partners of ai(1260) 
and 6i(1235), namely, K\a and Kib, respectively, are not mass eigenstates and they are mixed 
together due to the strange and non-strange light quark mass difference. We write 

Ki (1270) = K lA sin 6 Kl + K 1B cos 8 Kl , 

lfi(1400) = Kia cos 9k 1 — Kib sin 8k 1 ■ (2.1) 
If the mixing angle is 45° and {Kp\K\ B ) = (K p\K\a) , one can show that #i(1270) is allowed to 



decay into Kp but not K*ir, and vice versa for #i(1400) jl9l | 



From the experimental information on masses and the partial rates of iTi(1270) and #1(1400), 



Suzuki found two possible solutions with a two- fold ambiguity, {OrA ~ 33° and 57° [20(]. A similar 



constraint 35° < \6kxI ^ 55° is obtained in 2l| based solely on two parameters: the mass difference 
of the a% and b\ mesons and the ratio of the constituent quark masses. From the data of r — > 
#i(1270)^ T and #i(1400)^ T decays, the mixing angle is extracted to be ±37° and ±58° in [22]. As 
for the sign of the mixing angle, there is an argument favoring a negative 0k x ■ It has been pointed 
out in 23J] that the experimental measurement of the ratio of #17 production in B decays can be 



used to fix the sign of the mixing angle. Based on the covariant light-front quark model [24], it is 
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found 



m 1 

B(B ifi(1270)7) f 10.1 ± 6.2 (280 ± 200); for % = -58° (-37°) 



(2.2) 



B(B -> #i(1400)7) [ 0.02 ± 0.02 (0.05 ± 0.04); for 9 Kl = +58° (+37°). 

The Belle measurements B(B+ -> #+(1270)7) = (4.3 ± 0.9 ± 0.9) x 10~ 5 and B(B + -» 
#+(1400)7) < 1-5 x 10" 5 H clearly favor 9 Kl = -58° over 9 Kl = 58° and 9 Kl = -37° over 



9r 1 = 37°. In the ensuing discussions we will fix the sign of 9k x t° be negative. 

Likewise, the 3 Pi states /i(1285) and /i(1420) have mixing due to SU(3) breaking effects 

|/i(1285)) = \h)co S d aPl +|/ 8 }sin03 Pl , |/i(1420)) = -\f 1 )sin9 3Pl + |/ 8 )cos0 3Pl . (2.3) 

From the Gell-Mann-Okubo mass formula [l|, [3], it follows that 

4mt — rra 2 — 3rra 2 A , 19St .s 
cos 2 # 3Pl = — ^ , (2.4) 

3 l m /i(1420) _m /i(1285) 

where 

m K 1A = "^(1400) cos2 OKi + ^(1270) sin2 Oki • (2-5) 

Substituting this into Eq. ([23]) with Kl = -37° (-58°), we then obtain 0fpf = 27.9° (53.2°) and 
03p = 26.0° (52.1)° where the latter is obtained by replacing the meson mass squared m 2 by m 
throughout Eq. (|2.4p . The sign of the mixing angle can be determined from the mass relation [l| 

, — m 2 — 3m 2 f / nA9n \ 

tan03 Pl = Xl V ^ 5 /l(M20) . (2.6) 

2 V2(m 2 ai - m Ki J 

The previous phenomenological analyses suggest that 9z P ~ 50° [27J]. 2 Eliminating 9 from Eqs. 
(12.41) and (12.61) leads to the sum rule 



( m /!(i285) + m /i(i400))( 4m if 1A - m ai) - 3m5 l(1285) rnJ l(1400) = 8m KlA - 8m 2 KlA m 2 ai + 3m*. (2.7) 

This relation is satisfied for 3 P± octet mesons, but only approximately for l P\ states. Anyway, we 
shall use the mass relation (|2.4|) to fix the magnitude of the mixing angle and (|2.6|) to fix its sign. 

Since #*# and KKir are the dominant modes of /i(1420) whereas /o(1285) decays mainly to 
the 4-7T states, this suggests that the quark content is primarily ss for /i(1420) and ran for /i(1285). 
This may indicate that #3 P = 28° is slightly preferred. However, #3 Pl = 53° is equally acceptable. 



1 The sign of Qk x is intimately related to the relative sign of the Kia and K\b states. In the light-front 
quark model used in [23| and [24J, the decay constants of K\a and K\b are of opposite sign, while the 
B — > K\a and B — > #b form factors are of the same sign. It is other way around in the present work: 
the decay constants of K\a and K\b have the same signs, while the B — > K\a and £? — ► K\b form 
factors are opposite in sign. The two schemes are related via a redefinition of the K\a or K\b state, i.e. 
Kia —> —Kia or Kib —> —K\b- To write down Eq. (|2.2[) we have used our convention for Kia and K\b 
states. 

2 If a mixing angle 6a Pl of order 50° can be independently inferred from other processes, this will imply a 
preference of |0ki| = 58° over \0Ki \ — 37°. However, the phenomenological analysis in [27J is not robust 
and the Gell-Mann-Okubo mass formula employed there is not a correct one. 
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Similarly, for I 1 Pi states, /ii(1170) and /ii(1380) may be mixed in terms of the pure octet /ig 
and singlet h\, 

|fci(1170)) = |/ii)cos6»i Pl + \h 8 }sm6i Pl , \h x (1380)} = -|/n> sin 9x Pl + \h 8 ) cos 6i Pl . (2.8) 

Again from the Gell-Mann-Okubo mass formula, we obtain 

4ml- —ml —3m? 
cos ()i Pl = — -, — ^- , (2.9) 

3 I m h!(1380) ~~ m hi(1170) 



where 



m. 



l K 1B = m ^(i400) sin °K 1 + m Ki{1270) cos 6 Kl . (2.10) 

We obtain 6^ d = -18.1° (25.2°) and 0^ = 23.8° (-18.3°) for 6 Kl = -37° (-58°), where the sign 
of the mixing angle is determined from the mass relation 

4mi — m? —3m? d-,™-, 

, a J<1B 01 fti(1170) / 01 ,\ 

tan pi p, = — 7 =—z - — ^ - . (2.11) 



1Pl " 2y/2(ml-m 2 K 



ID ' 



B. Decay constants 

Decay constants of pseudoscalar and axial-vector mesons are defined as 

(P(p) 1 1 0) = -if P q» ( 3il) Pi(p, A)|g 2 7„759i|0) = if 3Pl{ i Pl) m3 Pli i Pl) e^>. (2.12) 
For axial-vector mesons, the transverse decay constant is defined via the tensor current by 

(V>I\(p, X)\q 2 a^ 5qi \0) = -f^^f - e -^) , (2.13) 

or 

< 3 «Pi(p, A)|g 2 ^?i|0) = -if k.)^ W e (A)/> ( 2 - 14 ) 

where we have applied the identity a a p^ = —^ a p l j, u a tlu with the sign convention eoi23 = 1- 
Since the tensor current is not conserved, the transverse decay constant / is scale dependent. 
Because of charge conjugation invariance, the decay constant of the l P\ non-strange neutral meson 
6^(1235) must be zero. In the isospin limit, the decay constant of the charged b\ vanishes due 
to the fact that the &i has even G-parity and that the relevant weak axial-vector current is odd 
under G transformation. Hence, f,± is very small in reality. Note that the matrix element of 
the pseudoscalar density vanishes, ( 3 ^Pi(p, £)|<?275<?i|0) = 0, which can be seen by applying the 
equation of motion. As for the strange axial vector mesons, the 3 Pi and 1 P± states transfer under 
charge conjunction as 

M b a ( 3 P l )^M^( 3 P 1 ), M b a ( 1 Pi)^-M^( 1 P 1 ), (a,b = 1,2,3). (2.15) 

Since the weak axial-vector current transfers as (A^) 1 ^ — > (Au)% under charge conjugation, it is clear 
that fi Pl = in the SU(3) limit [20]. By the same token, the decay constant f£p vanishes in the 
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SU(3) limit. Note for scalar mesons, their decay constants also vanish in the same limit, which can 
be easily seen by applying equations of motion to obtain 

m sfs = K m l ~ m 2)(0\qiq2\S), (2.16) 

with rrtj being the mass of the quark qi. 

The ai(1260) decay constant f ai = 238 ± 10 MeV obtained using the QCD sum rule method 



181 ] is similar to the p meson one, f p ~ 216 MeV. This means that the ai(1260) can be regarded as 



the scalar partner of the p, as it should be. To compute the decay constant fb 1 for the charged b\, 
one needs to specify the u and d quark mass difference in the model calculation. In the covariant 
light-front quark model 24], if we increase the constituent d quark mass by an amount of 5 ± 2 
MeV relative to the u quark one, we find = 0.6 ± 0.2 MeV which is highly suppressed. As we 
shall see below, the decay constant is related to the transverse one by the relation [see Eq. 

f bl =tiMal M (»), (2.17) 

where a{J' bl is the zeroth Gegenbauer moment of ^jj 1 to be defined later. The quantities and 
aQ bl can be calculated in the QCD sum rule approach with the results f^r = (180 ±8) MeV [3] (cf. 
Table U) and ajj' bl = 0.0028 ± 0.0026 for b{ at p = 1 GeV. (Note that for bf, al' bl has an opposite 
sign due to G-parity.) Again, is very small, of order 0.5 MeV, in agreement with the estimation 
based on the light-front quark model. In 14|, the decay constants of a\ and b\ are derived using the 
K 1A - Kib mixing angle 9 Kl and SU(3) symmetry: (/ 6lJ / 0l ) = (74, 215) MeV for Kl = 32° and 
(—28, 223) MeV for 6k 1 = 58°. It seems to us that the b\ decay constant derived in this manner is 
too big. 

Introducing the decay constants /| 1 ( 12 85) an d // 1 (i42o) 

(0|g 7 M75g|/i(1285)(P,A)) = -im /l(1285) /^ (1285) eW , (2.18) 
(0|S7mT5?|/i(1420)(P,A)) = ~i™ Ml4W ) fj^ofP > ( 2 - 19 ) 

we obtain 

/« = iK—Hth — CO s03 Pl + m /s s infe Pl =172 ±23 (178 ± 22) MeV, (2.20) 

// l( i285) = %-- ! ^cos03 Pl -%-^^sin03p 1 = -72±13(29±18)MeV, (2.21) 

V3 TOft (1285) V6 W/!(1285) 

// l( i420) = -% m/l sinfe Pl + % m/8 cos 6 3Pl = -55 ±10 (23 ±11) MeV, (2.22) 

V ; V 3^ (1420) V6 m /i(1420) 



//i(1420) 



l /l(1420) V m /i(1420) 

fh m h ^ a 2 ffs m h 



V / 3 m /i(1420) 1 \/6 "2/!(1420) 1 

= -219 ± 27 (-230 ± 26) MeV , (2.23) 

corresponding to 9s Pl = 53.2° (27.9°), where we have used the QCD sum rule results for and 
f fs |lj (see Table U). 

The decay constants for Ki(1270) and Ki(1400) defined by (with q = u or d) 

<0|g 7 M75s|Ki(1270)(P,A)> = -if K ,(i270)rn Kl{1270) eP, (2.24) 
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and 

(0|g 7 M75s|Ki(1400)(P,A)> = -i f Kl (uoo) m Kl[U00) eft (2.25) 
are related to fx 1A and fx 1B by 

/iCi(i270) = (fK 1A m KlA sin 9 Kl + fK lB m KlB cos 9 Kl ), 

m K 1 {1270) 

/jfifuoo) = " (fK 1A m KlA cos6 Kl - fK lB m KlB sva.9 Kl ). (2.26) 

m iCl(1400) 

Just as the previous b\ case, the decay constant fx 1B is related to the transverse one fx 1B by the 
relation /k 1b = fK 1B (^) a o' KlB 0-0 ■ ^ we apply the QCD sum rule results for /k- 1A) (see Table 
U) and aj'^ 113 (cf. Table |VJ) , we will obtain 

M(i270) = -137 ±15MeV, rV l(1400 ) = 199 ± 10 MeV, for Kl = -37°, 
M(i270) = -207 ±7 MeV, / Xl{1400) = 141 ± 14 MeV, for 9 Kl = -58°. (2.27) 

However, we would like to make two remarks. First, we do have the experimental information on the 
decay constant of i^i(1270). 3 From the measured branching ratio of r —* (l270)v T by ALEPH 



28], 0(r~ -> tff (1270)i/ T ) = (4.7 ± 1.1) x 1(T 3 , the decay constant of ifi(1270) is extracted to be 
221 



/^(mo) =175 ±19 MeV, (2.28) 



where use has been made of the formula 

g 2 F]v |2 ,2 K + Hi)W - m ^ 2 

1d7t m: 



r r - = ^|K,| 2 /i,— (2.29) 



Second, as pointed out in 2J|, the decay constants of 3 Pi have opposite signs to that of x Pi in 



the covariant light-front quark model. The large error with the QCD sum rule result of ajj' KlB = 
0.14 ± 0.15 is already an indication of possible large sum rule uncertainties in this quantity. 

In order to reduce the theoretical uncertainties with the K\ decay constant, we shall use the ex- 
perimental value of /r- 1 (i270) to fix the input parameters 0k 1a and 0k 1b appearing in the Gaussian- 



type wave function in the covariant quark model 241 ] . We obtain 



and 



0.375 GeV; for 9 Kl = -37°, 
Pk 1a = Pk 1b = \ Q 313 GeV . for % = _ 5g0) (2.30) 



f KlA = 293 MeV, f KlB = 15 MeV, for 9 Kl = -37°, 

/V 1A = 207 MeV, /^ ls = 12MeV, for 9 Kl = -58°. (2.31) 



3 The large experimental error with the -ftTi(1400) production in the r decay, namely, B(t~ — ► 
iff (1400)^) = (1.7 ± 2.6) x 10~ 3 [l|, does not provide sensible information for the ifi(1400) decay 
constant. 
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TABLE I: Summary of the decay constants fsp l and /ip (1 GeV) in units of MeV obtained from 
QCD sum rule methods [181 ] . 



3 PI 


ai(1260) 


h 


fs 


Kia 




238 ± 10 


245 ± 13 


239 ± 13 


250 ± 13 


'Pi 


&i(1235) 


hi 


/i 8 




Zip, 


180 ±8 


180 ± 12 


190 ± 10 


190 ± 10 



Therefore, we have 

M(1270) = -175 ± 11 MeV , / Xl (i400) = 235 ± 12 MeV , for 9 Kl = -37°, 

/^ l( i270) = "175 ± 15 MeV, /^(uoo) = 112 ± 12 MeV , for Kl = -58°. (2.32) 

In complete analogy to the discussion for l 3 Pi states, we introduce the tensor couplings for l 1 ^! 
states 

(0|^g|/ il (1170)(P,A)) = ifil q {im W e ?A)^, ( 2 -33) 

(0|^g|/ il (1380)(P,A)) = <f 1380 ) W e fA)^> ( 2 -34) 

and then obtain 

/ifiiTD) = % cos 01 Pi + % sin 01 ft = 75 ± 8 (127 ± 7) MeV , (2.35) 

/M1170) = % COS ^ ~ T^f Sin ^ = 147 ± 8 ( 28 ± 10 ) MeV ' ( 2 - 36 ) 

fhSsso) = -^sinftp 1 + ^cosft Pl = 106±5(26±7)MeV, (2.37) 

/M1380) = sinft A " ^=f cos 01 Pa = -115 ± 9 (-185 ± 10) MeV , (2.38) 

corresponding to flip = — 18.1° (25.2°) where we have used the QCD sum rule results for and 
ffc given in Table U Q • 

As for strange axial-vector mesons, we have (with q = u,d) 

(0|g < 7 AB/a |lfi(1270)(p,A)) = i/^ l(1270) W e fA)/ 

= sin6» Kl + /£ 1B cos fl^i) e^ef^p 13 (2.39) 

and 

<0|9<r^a|-Ki(1400)(p, A)> = ifk^uoo) ^p^V? 

= i{fk XA ™e Kl - fk 1B sm6 Kl ) e^ef A /. (2.40) 

As will be shown in Sec. II. D below, the decay constants Jk 1a anci fk 1A are related via 

fkM = fK 1A 4' KlA ^)- (2-4i) 



s 



From Tables U and [V] we obtain (at the scale /x = 1 GeV) 

/kx(1270) = 140 ± 22 MeV . /^(woo) = 130 ± 25 MeV , for % = -37°, 
/i(i270) = 8 4 ± 2 5 MeV , /i l(1400) = 172 ± 21 MeV , for 9 Kl = -58°. (2.42) 



Form factors 



The form factors for the B ^ A and B — > P transitions are defined as 

2 



(A(p,A)|A M |P(p B )) 

<A(p,A)|y M |5(pB)) 



-Um B - m A )$>V? A tf) - • pb)(pb + ^^^ 
I p m B — m A 

e (\)PB 



(P(p)\V,\B( PB )) 



(PB + P) A 



BA(J2\ 



m B — mj 



Ff P (<?) + 



' v BP (n 2\ 



(2.43) 



(2.44) 



where q = p B - p, V 3 BA (0) = V BA (0), Ff F {0) = F BF {0) and 

2myi 2myi 
In the literature the decay constant and the form factors of the axial vector mesons are often 
defined in different manner. For example, in 24] they are defined as 4 



<A(p,A)|A M |0> = f A m A e 



(A)* 



(2.45) 



(A(p,X)\A^\B(p B )) 
{A{p,\)\V^\B{p B )) 



m B - rriA 
-i\ (m B 



e, ua ^ x) p a B p^A BA (q 2 ), 

mA^V^tf) - ^UpbKpb+p),^^ 



(2.46) 



It has been checked in the covariant light-front quark model that the form factors V 1 2 1 {q 2 ) and 
A B3pi (q 2 ) defined in Eq. (j2.46j) are indeed positively defined. We would like to ask if the B — > 3 Pi 
transition form factors defined in Eq. (|2.43p are also positively defined. This can be checked by 
considering the factorizable amplitudes for the decay B — > AP 



x (ba,p) = {Pmv _ A)^\0)(A(p)\(V - Ar\B( PB )), 
X (BP,A) = ( A (q)\(V - AU0)(P(p)\(V - Ar\B( PB )). 



(2.47) 



4 Since the convention g 1234 = +1 is adopted in |24( while £1234 = +1 is used in the present work, we have 
put an additional minus sign for the matrix element (A(p, A)|A M |S(ps)). 
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TABLE II: Form factors for B tt, K, ai(1260), 61 (1235), K 1A , K 1B 



covariant light-front model 2J| are fitted to the 3-parameter form Eq 
factor V2 denoted by * for which the fit formula Eq. (|2.5ip is used. 



transitions obtained in the 
(|2,50|) except for the form 



F 


F(0) 


^(9max) 


a 


b 


F 


F(0) 




a 


b 


F m 


0.25 


1.16 


1.73 


0.95 


Ft 


0.25 


0.86 


0.84 


0.10 


pBK 


0.35 


2.17 


1.58 


0.68 


F BK 


0.35 


0.80 


0.71 


0.04 




0.25 


0.76 


1.51 


0.64 


-trBai 
v 


0.13 


0.32 


1.71 


1.23 


yBai 


0.37 


0.42 


0.29 


0.14 


v 2 Bai 


0.18 


0.36 


1.14 


0.49 


A B bl 


-0.10 


-0.23 


1.92 


1.62 




-0.39 


-0.98 


1.41 


0.66 


yB bl 


0.18 


0.36 


1.03 


0.32 


V B * 


0.03* 


-0.15* 


2.13* 


2.39* 




0.26 


0.69 


1.47 


0.59 




0.14 


0.31 


1.62 


1.14 


yBK 1A 


0.39 


0.42 


0.21 


0.16 


yBK 1A 


0.17 


0.30 


1.02 


0.45 


ABK 1B 


-0.11 


-0.25 


1.88 


1.53 


\rBK\s 


-0.41 


-0.99 


1.40 


0.64 


yBK 1B 


-0.19 


-0.35 


0.96 


0.30 


yBK 1B 


0.05* 


0.16* 


1.78* 


2.12* 



TABLE III: Form factors for B -» oi(1260), h (1235), K 1A , K 1B , fu / 8 , hi, h 8 transitions at q 2 = 
obtained in the framework of the light-cone sum rule approach 31[. Uncertainties arise from the 
Borel window and the input parameters. 



v Bai (o) 


n ono+0.022 
<J.OU<5_o.o35 


v Bbl (o) 


^ n +0 039 
u.ooo_ 003 3 


V BKlA (0) 


S16+ - 048 
u.oio_ 042 


V BKxB (0) 


-0 360+ 030 

U.ODU_q 02g 


v Bh (o) 


n 1 o-| +0.018 
U - 18i -0.021 


v Bhl (o) 


-0 214+ 021 

U - Zi ^-0.012 


V Bf »(0) 


n i24+ - 015 
u.iz4_ 004 


V Bh *(0) 


n 1 ^S+0-0!6 



We obtain 



X (BA,P) = -2if P m A V BA (q 2 )(e* {x)PB ) , X^ BP ^ = -2if A m A F BP (q 2 )(e* w p B ) , (2.48) 

from Eqs. (12U21) and fl!Q3l) and 

= 2f P m A V BA (q 2 )(e* w p B ) , = -2f A m A F BP V)(<a)Pb) , (2.49) 

from Eqs. (|2.45p and (|2.46p . Since f A for the 3 Pi meson is negative in the light-front model 
calculation (see Eq. (2.23) and Table III in 2Jj), the relative sign between X^ PA ^ and X^ P ' A ^ 
is positive. This means that the relative sign in Eq. (12.480 is also positive provided that the decay 
constant f A and the form factor Vn A defined in Eqs. (I2.12P and (12.43D . respectively, are of the 
same sign. Indeed, it is found in [16j that if f A is chosen to be positive for the ai(1260) meson, the 
form factor V Bai is indeed positive according to the sum rule calculation. 

The form factors for B — > tt, K, Gq(1260), b\ (1235), K\ A , K\ B transitions have been calculated 
in the relativistic covariant light-front (CLF) quark model 24| (Table [II]) 5 and in the framework 



As explained in the footnote before Eq. (|2.2[) . we need to put additional minus signs to the B — > 1 P 1 form 
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of the light-cone sum rule (LCSR) approach [31]. In the CLF model, the momentum dependence 
of the physical form factors is determined by first fitting the form factors obtained in the spacelike 
region to a 3-parameter function in q 2 and then analytically continuing them to the timelike region. 
Some of the V2(q 2 ) form factors in P —> A transitions are fitted to a different 3-parameter form so 
that the fit parameters are stable within the chosen q 2 range. 

Except for the form factor V2 to be discussed below, it is found in [241 ] that the momentum 
dependence of form factors in the spacelike region can be well parameterized and reproduced in the 
three-parameter form: 

F{q2) = l-a(q 2 /m 2 B )+b(q 2 /m 2 B ) 2 ' (2 - 50) 

for B — > M transitions. The parameters a, b and F(0) are first determined in the spacelike region. 
We then employ this parametrization to determine the physical form factors at q 2 > 0. In practice, 
these parameters are generally insensitive to the q 2 range to be fitted except for the form factor 
V2(q 2 ) in B — > 1 P\ transitions. The corresponding parameters a and b are rather sensitive to the 
chosen range for q 2 . This sensitivity is attributed to the fact that the form factor Vz(q 2 ) approaches 
to zero at very large — \q 2 \ where the three-parameter parametrization (|2.50p becomes questionable. 
To overcome this difficulty, we will fit this form factor to the form 

F (J\ = EM ( 2 51) 

{q ' (l-q 2 /m 2 B )[l-a(q 2 /m 2 B )+b(q 2 /m 2 B ) 2 ] V " ' 

and achieve a substantial improvement 241 ]. 

Momentum dependence of the form factors calculated using the LCSR method is not shown in 
Table IIIH Since the pseudoscalar mesons considered in the present work are the light pion and the 
kaon, the form-factor q 2 dependence can be neglected for our purposes. 



TABLE IV: B — ► ai(1260) transition form factor V Bai at q 2 = in various models, where the QSR1 
is the traditional QCD sum rule approach and the QSR2 is the light-cone sum rule approach. 





CLF [24] 


ISGW2 [30] 


CQM [29] 


QSR1 [32] 


QSR2 [31] 


v Bai (o) 


0.13 


1.01 


1.20 


0.23 ± 0.05 


n oqo+0.022 
U.OUO_ 035 



In principle, the experimental measurements of ~W — * a~^ir~ and B° — > 6^7r~ will enable us to 
test the form factors V Bai and V$ bl respectively. There are several existing model calculations 
:br B — ► a\ form factors: one in a quark-meson model (CQM) 29J, one in the ISGW2 model 



30j, one in the light-front quark model [24j and two based on the QCD sum rule (QSR) [3l|, [32]. 
Predictions in various models are summarized in Table ITvl and in general they are quite different. 
For example, Vo Bai (0) obtained in the quark-meson model, 1.20, is larger than the value of the 



sum-rule prediction in [32|, 0.23 ± 0.05. If ai(1260) behaves as the scalar partner of the p meson, 



factors in Table ITT1 since in the convention of the present work, the form factors Pl and V® 

different signs. 



have 
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V Bai is expected to be similar to Aq P , which is of order 0.28 at q 2 = 24]. Indeed, the sum rule 
calculation by one of us (K.C.Y.) yields V Bai = 0.303±g;g|§ . Therefore, it appears to us that a 
magnitude of order unity for V Bai (0) as predicted by the ISGW2 model and CQM is very unlikely. 



The BaBar measurement of B u — ► cl^tt [4j] favors a value of V Bai (0) ~ 0.30, which is very close 
the LCSR result shown in Table [TTTJ 

Various B — > A form factors also have been calculated in the Isgur-Scora-Grinstein-Wise ( 
model 



SGW) 



24j, 



m 



3d . |33J based on the nonrelativistic constituent quark picture. As pointed out in 
general, the form factors at small q 2 in CLF and ISGW models agree within 40%. However 



BD* BD 1 ^ A ' 

Fq (q 2 ) and V 1 1 (q 2 ) have a very different q 2 behavior in these two models as q 2 increases. 
Relativistic effects are mild in B — > D** transitions but can manifest in heavy-to-light transitions 
at maximum recoil. For example, V^ ai (0) is found to be 0.13 in the CLF model, while it is as big 
as 1.01 in the ISGW2 model. 



.1/2 



D. Light-cone distribution amplitudes 



For an axial-vector meson, the chiral-even LCDAs are given by 



(A(P, \)\qi(y) W q 2 (x)\0) = im A £ due^ + ^ ^ *„(«) + e^* gf '(«)} , 

(2.52) 

(A(P, A)|gi(i/) 7M g2(x)|0) = -im A W ^p"** J* due^ +a ^ 2±^-, (2.53) 
with u (u = 1 — u) being the momentum fraction carried by q\{q~2), and the chiral-odd LCDAs read 

(A(P, X)\qi(y)a^ 75 q 2 (x)\0) = [' due^ + ^ \{e%*p v - e^* Pfl ) * ± (u) 



o 



+ (Pm^ - Pv*n) h f ( u ) I ( 2 - 54 ) 



(A(P, X)\qi(yh 5 q2(x)\Q) = m 2 A e^> z / due< u ™ + ^ -L— . (2.55) 



o 



Here, throughout the present discussion, we define z = y — x with z 2 = and introduce the 
light-like vector p^ = P^ — m\z^/(2Pz) with the meson's momentum P 2 = m\. Moreover, the 
meson polarization vector e M has been decomposed into longitudinal (cjj^ ) an d transverse (ejj^*) 
projections defined as 



e ii M - Pz ^ p z ^ r ±>m ~ M h** ' 1 j 



respectively. The LCDAs <J>||,3>j_ are of twist-2, and , hfc of twist-3. Due to G-parity, 

3?[|jSj^ anc ^ #± are s y mme t r ic (antisymmetric) with the replacement of u — > 1 — u for 3 Pi ( 1 i 3 i) 
states , whereas 3>j_, /in and /i^ are antisymmetric (symmetric) in the SU(3) limit [la ]. We restrict 
ourselves to two-parton LCDAs with twist-3 accuracy. 
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Assuming that the axial-vector meson moves along the z— axis, the derivation for the light-cone 
projection operator of an axial- vector meson in the momentum space is in complete analogy to 
the case of the vector meson. We separate the longitudinal and transverse parts for the projection 
operator: 

M s A a = M s A all +M s A a± , (2.57) 
where only the longitudinal part is relevant in the present study and is given by 



Mi 



.1 m A{e* (X ) n +) , , s im A m A(e* (X ) n +) 

-l — <P||(u) — 

4 2 ' lb ^ ' 4 2E 



^v^lbn-n+h^ {u) 



+ 



iE ( U dv ($,» - h\!\v)) a^nt 
Jo 11 



/if>)(t 



dk± u 



75" 



k=up 



with the momentum of the quark q\ in the A meson being 



= uEnt + k^ + 



kA 



AuE 



(2.58) 



(2.59) 



for which E is the energy of the axial-vector meson and the term proportional to k\ is negligible. 



Here, for simplicity, we introduce two light-like vectors = (1,0,0,-1) and 



[1, 0,0,1). In 



general, the QCD factorization amplitudes can be recast to the form Jq duTi(M A • • •). 
The LCDAs <S> A ± (u) can be expanded in terms of Gegenbauer polynomials of the form: 



a Jj(^ + Ea ll(^f/ 2 (2.-l) 



(2.60) 



where the relevant decay constants in the above equation will be specified later. In the following 
we will discuss the LCDAs of 1 P\ and 3 Pi states separately: 



1. P\ mesons 

For the Q^AJu), due to the G-parity, only terms with odd (even) Gegenbauer moments survive 

in the SU(3) limit. Hence, the normalization condition for the twist-2 LCDA can be chosen 
as 

j\u^(u) = f 1 A Pi . (2.61) 
In the present work, we consider the approximation 

^(u) = f^&uu jl + 3a^' 1Pl (2n - 1) + 4' 1Pl | 
Likewise, we take 

$J Pl (u) = / 1 - i Pi 6nn|a||' 1Pl +3a|' 1Pl (2n-l)+o5' 1Pl ^ 
with the normalization condition 

j\u^(u)=fi Pl . (2.64) 



5(2u - l) 2 - 1 



(2.62) 



5(2u- l) 2 - 1 



(2.63) 
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This normalization together with Eq. (|2.52p leads to Eq. (|2.12p for the definition of the 1 P± meson 
decay constant. Eqs. (|2.63|) and (|2,64|) lead to the relation 

fi Pl =fi A Pl ^)al' 1Pl (n). (2.65) 

i II 1 Pi 

The scale dependence of /fp must be compensated by that of the Gegenbauer moment a ' to 

ensure the scale independence of /ip 1 . In principle, we can also use the decay constant fi Pl to 

construct the LCDA ^m^ 1 - However, fi Pl vanishes for the neutral 61 (1235) and is very small for 

the charged &i(1235). This implies a vanishing or very small &\\ Pl unless the Gegenbauer moments 

I 1 -Pi 

cii are very large. Hence, it is more convenient to employ the non- vanishing decay constant 
f^ Pi to construct ^y^ 1 - This is very similar to the scalar meson case where the twist-2 light-cone 
distribution amplitude <&s is expressed in the form 



$s(x,/x) = f s (jJk)6x(l-x) 



m=l 
'Pi 



with fs being defined as (S\q2qi\0) = msfs- Now <&,, 1 can be recast to the form 



^(u) = / 1 p 1 6^|l + / i 1 p 1 ^af' 1Pl cf 2 (2n-l)|, 



(2.66) 



(2.67) 



with fj,i Pl = l/a{j' 1 . For the neutral 61 (1235), fb 1 vanishes and fj,^ becomes divergent, but the 
combination fbi^bi ls finite [35:]- Recall that for the scalar meson case, its LCDA also can be 
expressed in the form 



^s(x,n) = fs6x(l - x) 



1 + MS 



m=l 



(2.68) 



where fs = fJ-sfs and the equation of motion leads to fig = ms/(ffl2(fi) — mi (//)). However, unlike 
the case for scalar mesons, the decay constants /fp and fi Pl cannot be related by equations of 
motion. 

When the three-parton distributions and terms proportional to the light quark masses are ne- 
glected, the twist-3 distribution amplitudes can be related to the twist-2 3>j_(ii) for the transversely 
polarized axial- vector meson by Wandzura-Wilczek relations fl~8| : 



hf{v) = {2v-l) 



du 
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1 $±{u) 



du 



h'PHv) 



v <$>±(u) 



du(<5>±(u) - hf\u)) = vv 



du 







1 $ ± («) 
u 

v <S>±{u) 
u 



u 

■du 
du 



(2i; -!)*«(«) , 



-2* a (t;) 
1 $ ± (u) 



du 



u 



VV$ a (v) . 



The twist-3 LCDA <3? a (tt, n) satisfies the normalization 

/ $ a (u)du = 0, 
Jo 

and has the general expression 



(2u-l) + f^ay p i(riP n+1 (2u-l) 



n=i 



(2.69) 

(2.70) 
(2.71) 



where P n {u) are the Legendre polynomials. 
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2. 3 Pi mesons 



In analogue to the 1 P± case, we consider the approximations: 



$ Pl (w) = f 3Pl 6uu( 1 + 3af " 1 (2n-l)+a 



ft, 



5(2u- l) 2 - 1 



Sjf^u) = / 3Pl 6mt <{ Oq rjri + 3a] Lr ' ri (2n - 1) + 



-L, -Pi 



3 A 3 



X, Pi 



5(2« - l) 2 



(2.72) 
(2.73) 



In the SU(3) limit, only terms with even (odd) Gegenbauer moments for ^iifn survive due to the 



1 1 3 p , i , ^ ^ 

G-parity. Hence, and a 2 1 vanish in the SU(3) limit. The LCDAs respect the normalization 

conditions 



-, 3 Pi 



1 1 du$! Pl (u) = fs Pl , ( l du^ (u) = f 3 \ , 
Jo 11 JO 



and 



duhf{u) = J duh^'(u) = 0. 
The latter is valid in the SU(3) limit. Therefore, we obtain 

/spiCaO = /3PiOo' 3Pl (Ai), 

and 

*^(«) = / 3 - L Pi 6nn{l + ^p 1 ^a^ 3pi Cf /2 (2n-l) 



Mi 



i=l 



with p3p 1 = l/a ' 1 - The twist-3 LCDA <3? a has the expression 



$ a Pl (u, Ai )=3/ 3Pl 



z ± ' 3Pl (2u - 1) + E 4' 3 *(/i)P n+1 (2u - 1) 



n=l 



(2.74) 
(2.75) 

(2.76) 

(2.77) 

(2.78) 



Most of the relevant Gegenbauer moments a. 



method 



have been evaluated using the QCD sum rule 



181 ] . The results are summarized in Table El 



For the pseudoscalar meson LCDAs we use 

$ P (tt) = fpQuul l + 3af(2it- 1) +a P - 



5(2u - l) 2 - 1 



= f P , 



6 



/p«(1 - it), 



(2.79) 



where <E> P and $0- are twist-3 LCDAs. We shall employ the sum rule results for the Gegenbauer 
moments of pseudoscalar mesons 

/i = 1.0GeV: of = 0.06 ± 0.03 , of = 0.25 ± 0.15, = 0, of = 0.25 ± 0.15 , 
/ u = 2.1GeV: of = 0.05 ± 0.02 , of = 0.17 ± 0.10, a\ = 0, a£ = 0.17 ± 0.10. (2.80) 
Note that in this paper the G-parity violating parameters (of , a\' KlA , %'^ 1A , a^' KlB and o)q2 1B ) 
are for mesons containing a strange quark. For mesons involving an anti-strange quark, the signs of 
G-parity violating parameters have to be flipped due to the G-parity. The integral of the B meson 
wave function is parameterized as 37J 



' ^ *f (P) = ^ 



(2.81) 



/o 1 - P A B 

where 1 — p is the momentum fraction carried by the light spectator quark in the B meson. Here 
we use A B (1 GeV) = (350 ± 100) MeV. 



15 



TABLE V: Gegenbauer moments of and $n for l 3 Pi and \ l P\ mesons, respectively, taken from 



18] 





||, ai (1260) 

a 2 


11^ 

«2 1 


IU 8 3P1 

«2 8 




4*** 


1 GeV 


-0.02 ±0.02 


-0.04 ± 0.03 


-0.07 ±0.04 


-0.05 ± 0.03 


0.00 ±0.26 


2.2 GeV 


-0.01 ± 0.01 


-0.03 ± 0.02 


-0.05 ± 0.03 


-0.04 ± 0.02 


0.00 ± 0.22 




±,01(1260) 




V' 8 


-1 Kl A 


±,Kia 
V 




1 GeV 


-1.04 ±0.34 


-1.06 ±0.36 


-1.11 ±0.31 


-1.08 ±0.48 


0.08 ± 0.09 


0.02 ±0.20 


2.2 GeV 


-0.83 ± 0.27 


-0.84 ± 0.29 


-0.90 ± 0.25 


-0.88 ± 0.39 


0.07 ±0.08 


0.01 ±0.15 




II, bi (1235) 
"1 


ii^: pi 


af 8 


ll.-Kl R 


J,K 1B 




1 GeV 


-1.95 ±0.35 


-2.00 ± 0.35 


-1.95 ±0.35 


-1.95 ±0.45 


0.14 ±0.15 


0.02 ±0.10 


2.2 GeV 


-1.61 ±0.29 


-1.65 ±0.29 


-1.61 ±0.29 


-1.57 ±0.37 


0.14 ±0.15 


0.01 ±0.07 




±,6i(1235) 

a 2 


±,h] Pl 
a 2 


±,hl Pl 
a 2 8 






1 GeV 


0.03 ±0.19 


0.18 ±0.22 


0.14 ±0.22 


-0.02 ±0.22 


0.17 ±0.22 


2.2 GeV 


0.02 ±0.15 


0.14 ±0.17 


0.11 ±0.17 


-0.02 ±0.17 


0.14 ±0.18 



E. Other input parameters 



For the CKM matrix elements, we use the Wolfenstein parameters A = 0.818, A = 0.22568, 
p = 0.141 and f\ = 0.348 38]. The corresponding three unitarity angles are a = 90.0°, (3 = 22.1° 
and 7 = 68.0°. 

For the running quark masses we shall use 

m b {m h ) = 4.2 GeV, m 6 (2.1 GeV) = 4.95 GeV, m 6 (l GeV) = 6.89 GeV, 

m c (m b ) = 1.3 GeV, m c (2.1 GeV) = 1.51 GeV, 

m s (2.1GeV) = 90MeV, m s (lGeV) = 119 MeV, 

m d (lGeV) = 6.3 MeV, m u (lGeV) = 3.5 MeV. (2.82) 
The uncertainty of the strange quark mass is assigned to be m s (2.1 GeV) = 90 ± 20 MeV. 
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III. B — ► AP DECAYS IN QCD FACTORIZATION 



We shall use the QCD factorization approach [37J, 139J to study the short-distance contributions 
to the B AP decays with A = oi(1260), /i(1285), /i(1420), Id (1270), 61 (1235), hi (1170), 
hi(1380), Ki(1400) and P = ir,K. It should be stressed that in order to define the LCDAs of 
axial-vector mesons properly, it is necessary to include the decay constants. However, for practical 
calculations, it is more convenient to factor out the decay constants in the LCDAs and put them 
back in the appropriate places. Recall that ^y^ 1 has two equivalent expressions, namely, Eqs. (|2.63p 

and (|2.67p . However, we found out that it is most convenient to use Eq. (|2.63p for the LCDA &\\ Pl 
which amounts to treating the axial-vector decay constant of 1 P\ as fi Pl - (Of course, this does not 
mean that fi Pl is equal to fi~ Pl -) Likewise, we shall use Eq. (|2.77p rather than Eq. (|2.73p for the 
LCDA 

In QCD factorization, the factorizable amplitudes of above-mentioned decays are collected in 
Appendix A. They are expressed in terms of the flavor operators a? and the annihilation operators 
6f with p = u,c which can be calculated in the QCD factorization approach [37j|. The flavor 
operators aP i are basically the Wilson coefficients in conjunction with short-distance nonfactorizable 
corrections such as vertex corrections and hard spectator interactions. In general, they have the 
expressions 



37 



39] 



af(MiM 2 ) 



Ci + 



Ci±\ 



+ 



N r 4tt 



Ni(M 2 ) 
Vi(M 2 



<$>f 2 (x)dx 



+ 



4vr^ 



H i (M 1 M 2 ) +Pf(M 2 ), 



(3.1) 



where i = 1, • • • , 10, the upper (lower) signs apply when i is odd (even), c, are the Wilson coefficients, 
C F = (iV c 2 - 1)/(2JV C ) with N c = 3, M 2 is the emitted meson and Mi shares the same spectator 
quark with the B meson. The quantities Vi{M 2 ) account for vertex corrections, Hi(M\M 2 ) for hard 
spectator interactions with a hard gluon exchange between the emitted meson and the spectator 



quark of the B meson and Pi{M 2 ) 
Ni(M 2 ) reads 



for penguin contractions. The expression of the quantities 



Ni(M 2 



i = 6, 8 and M 2 
else. 



A, 



(3.2) 



0, 

Note that Ni(M 2 ) vanishes for i = 6, 8 and M 2 = A as a consequence of Eq. (|2.70p . The subscript 
|| of $ in the first term of Eq. (|3.ip reminds us of the fact that it is the longitudinal component of 
the axial-vector meson's LCDA that contributes to the B — > AP decay amplitude. Specifically, we 
have 



du$l Pl {u) 



1. 



(3.3) 



where we have factored out the decay constant (f3 Pl ) of (<I> 
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The vertex corrections in Eq. (J3JJ) are given by 

" 1 m b 



Vi(M 2 ) = { 



dx$M 2 (x) 12 In — - 18 + g{x) 
p 

dx®M 2 [-121b— + 6-g(l-x) 
L p 



dx<& rri2 (x) — 6 + h(x) 



(i = 1-4,9, 10), 
(t = 5,7), 
(i = 6,8), 



(3.4) 



with 



9{x) 



3 ( — In x — m 



+ 



1 — x 



2Li 2 (x) -ln 2 x + 



2lnx 
1 — x 



(3 + 2z-k) In x — (x <-> 1 — x) 



h(x) = 2Li 2 (x) — In 2 x — (1 + 2«7r) lnx — (x «-> 1 — x) 



(3.5) 



where 3>m (^m) is the twist-2 (twist-3) light-cone distribution amplitude of the meson M. More 
specifically, $m = <&p, <f> m = <£ p for M = P and $m = $ m = $ s for M = i. For the general 
LCDAs 

oo 

a + Y,an{lj)Cl /2 {2x-l) , 



?1=1 



$m0*0 = 6x(l - x) 

oo 

$m(x) = P + Y,M^)Pn(2x-l), 



(3.6) 



n=l 



the vertex corrections read 



< . / , m 6 1 \ /ll \ 21 /79 2ivr\ , , 

Vi(M) = \12 In -i - 18 - - - 3ivrj « + - 3ivrj ai - — a 2 + ^— - —J a 3 + • • • , (3.7) 

for i = 1-4,9,10, 

, , f rn b 1 \ /ll \ 21 /79 2ivr\ , , 

Fi(M) = (-12111 — + 6 + - + 3zttJ ao + (^y - 3zttJ ttl + — a 2 + + —J a 3 + • • • , (3.8) 



for i = 5,7, 



Vi(M) = — 6/3 + (3 — i2n)(3 1 + 



19 ivr 
18 ~ Y 



/?3 + 



(3.9) 



for i = 6,8. 

As for the hard spectator function H , it has the expression 



fli(MiM 2 ) 



X (BAh,M 2 ) J p 



-21$ B (p) / ^$ A/2 (£) / _L 



o £ 



1 cfy 



o ?? 



<Kw,(//)±rf | <&„„(>/> 



(3.10) 



for i = 1 — 4, 9, 10, where the upper sign is for Mi = P and the lower sign for M\ = A, 



X(BM 1 ,M2) Jo p JO s JO ^ 



$ Ml (7?)±r^|$ mi (r7) 



(3.11) 
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for i = 5, 7 and Hi = for i = 6,8, £ = 1 — £ and rj = 1 — rj. In above equations, 



2mp 



2771^4 



x m j ,(/i)(m 2 + rai)(/i)' x m 6 (/i) ' 

and x^ BMl ' M2 ^ is the factorizable amplitude defined in Eq. (12.4Tf) . 

Weak annihilation contributions are described by the terms 6j, and &i,EW hi Eq. 
have the expressions 



(3.12) 



11) which 



h 

&3,EW 



c x A\, 



Cf 
Cf 



Cf r 

iVf 
Cp 
AT 2 



c 4 A + c 6 ^^ 



63 

c 9 A + c 7 (A| + A{) + 7V C C 8 A 
cio4 + c 8 A, 



c 3 A[ + c 5 (Ai + A{)+N c c 6 A f 3 



(3.13) 



where the subscripts 1,2,3 of A^f denote the annihilation amplitudes induced from (V — A)(V — A), 
(V — A)(V + A) and (S — P)(S + P) operators, respectively, and the superscripts i and / refer to 
gluon emission from the initial and final-state quarks, respectively. Their explicit expressions are: 



A\ 



A\ 



A{ 



$p(x)$ A (y) 
<f> A (x)$ P (y) 

$p{x)$ A (y) 

$ A (x)$p(y) 

-r£<S>p(x)<S> a (y) 
r P <f> A (x)%(y) 



y(l-xy) 


+ 


Wy 


1 


+ 


1 


y(l-xy) 


Wy 


1 


+ 


1 


x(l-xy) 


ly 7 


1 


+ 


1 



r x r x $ P (x)$a(y) h 



xy(l-xy) 



A P 

"x T X 

r^r P ^ a (x)^> p {y) 

r*r P %{x)<5> a (y) £ 
r£r?$ a (x)%(y) ± 

r P %(x)^ A (y) Iyl j^ 



xy(l-xy) 



2.T 



xy(l-xy) 



-r$* P (x)* a (y) 2 -^f + r P %{x)$ A {y) HilM 



xy 



r P ® A {x)%{y) ^ - ti* a {x)* P {y) ^ 



for MiM 2 = AP, 
for M X M 2 = PA, 

for MxM 2 = AP, 
for M 1 M 2 = PA, 

for MiM 2 = AP, 
for MiM 2 = PA 

for MiM 2 = AP, 
for MiM 2 = PA, 



A[ = A{ = 0, 



(3.14) 



where [ • • • = ira s Jq 1 dxdy, x = 1 — x and y = 1 — y. Note that we have adopted the same convention 
as in [39(] that M\ contains an antiquark from the weak vertex with longitudinal fraction y, while 
M2 contains a quark from the weak vertex with momentum fraction x. 

Using the asymptotic distribution amplitudes for <£p, <3? p , ^^ Pl , ^ a Pl and the leading contribu- 
tions to $u Pl , 



<&p(u) = 6uu, 

= 1, < Pl (n) 
we obtain from Eq. (|3.14j) that 



(u) = 6uu, 



3a^ Pl (6^-6tt + l), $o^(u) = 3(2n - 1 



18a|' 1Pl ii'u(27j - 1), 
'Pi, 



(3.15) 



A/ 3 



PiP) 



3 X A - 4 + 



7T 



^ 3 %^rfX A (X A -3) 
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67ra s 
67ra s 
67ra s 
67ra s 

67ra s 



3 [X A - 4 + 



4( 3 PiP) 
4('PiP) 

4Cap) 

4( 3 PiP) « 6vra s (2X A -l) 



-3al' 1Pl (X A + 29 - 3tt 2 ) + r^r^MXi 

2- 



7T 



-3a\' 1Pl (3X A + 4-tt 2 )+ r^r p X A {X A - 2) 



Y 2 



2X A + 



3a 





2X A 



2X A + 4 - — 



4(^1 P) 



t S' 1Pl <X A (6X A - 11) + 2,r 1 x ^{2X A - 1){X A - 2) 



(3.16) 



and 



^l(P 3 Pi) = A\( 3 P 1 P), 
,4 i 2 (P 3 P 1 ) = 4( 3 P 1 P), 



M(p 1 Pi) 



-aU 1 



PiP), 



-AU l PiP), 



Al(p- S p 1 ) = -Ai( s p 1 p), 



A f 3 (P*P l ) = AiCP 1 P) 



ff3, 



2K A 1 U — -^H 

A\{P 1 P 1 )=A\{ 1 P 1 P), 

4(p i p 1 ) = -4( i p 1 p), 



(3.17) 



where the logarithmic divergences occurred in weak annihilation are described by the variable X A 

(3.18) 



' x A , 

o u 



1 Inn 1 ^ 
o u 2 



Following 37], these variables are parameterized as 

X A = ln 



(3.19) 



with the unknown real parameters p A and (fi A . Likewise, the endpoint divergence Xh in the hard 
spectator contributions can be parameterized in a similar manner. Following 34|, [40|, we adopt 
PA,H < 0.5 and arbitrary strong phases 4> Ai h with p At H = by default. 

Besides the penguin and annihilation contributions formally of order 1/nib, there may exist 
other power corrections which unfortunately cannot be studied in a systematical way as they are 
nonperturbative in nature. The so-called "charming penguin" contribution is one of the long- 
distance effects that have been widely discussed. The importance of this nonpertrubative effect 
has also been conjectured to be justified in the context of soft-collinear effective theory 41[. More 
recently, it has been shown that such an effect can be incorporated in final-state interactions j^ . 
However, in order to see the relevance of the charming penguin effect to B decays into scalar 
resonances, we need to await more data with better accuracy. 
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IV. NUMERICAL RESULTS 



A. Branching ratios 

The calculated branching ratios for the decays B — * An, AK with A = 
ai (1260), 61 (1235), ^(1270),^! (1400) ,/i (1285), /i (1420), /ti (1170), /ii (1380) are collected in Ta- 
bles IVmVIHl For B — > A transition form factors we use those obtained by the sum rule approach, 
Table [ffll The theoretical errors correspond to the uncertainties due to variation of (i) the Gegen- 
bauer moments (Table [Vj, the axial- vector meson decay constants, (ii) the heavy-to-light form 
factors and the strange quark mass, and (iii) the wave function of the B meson characterized by 
the parameter Xb, the power corrections due to weak annihilation and hard spectator interactions 
described by the parameters pa,h, 4>a,h, respectively. To obtain the errors shown in Tables IVII - 
IVIIIl we first scan randomly the points in the allowed ranges of the above seven parameters in 
three separated groups: the first two, the second two and the last three, and then add errors in 
each group in quadrature. 



1. B — > aiTT, a±K decays 

From Table IVT1 we see that the predictions for B° — > afir^ are in excellent agreement with the 
average of the BaBar and Belle measurements 0, 0]. BaBar has also measured time-dependent 
CP asymmetries in the decays 5° — > af-x^ {-djj. Using the measured parameter AC (see Sec. IV. B), 
BaBar is able to determine the rates of B° — ► a^ir~ and B° —* a^ir + separately, as shown in Table 
IVII It is expected that the latter governed by the decay constant of ai(1260) has a rate larger 
than the former as f ai S> f n . Again, theory is consistent the data within errors. However, there 
are some discrepancies between theory and experiment for a^vr - and cl^tt® modes. It appears that 
the relations B(B~ -» a?7r - ) > B(B° afir-) and B(B~ a^vr ) > B(B° -> a{ir + ) observed 
by BaBar are opposite to the naive expectation that B(B~ — ► a^Tr - ) < B(B° — > a^7r~) and 
B(B~ — ► a^7r°) < B(B° — > ajfvr + ). As for B — ► a\K decays, although the agreement with the 
data for 5° — > a\K~ is excellent, the expectation of B(B~ — ► a^X°) ~ B(B° — » a^i^ - ) is not 
consistent with experiment. More specifically, in QCDF we obtain the ratios 

_ g(g -» a°7r ) S o+0.02+0.06+0.12 / , . 1 (i7 , n7s \ 

R\ = — =0 7 — : - u -&J_ .o2-o.05-o.i7 (expt . 1.67±U.78), 

D(i3 — > aj"vr ) 

R * S !l^o - ^ = 0.62±g^±g^±8:Sg (expt : 1.26 ± 0.46), 

i? 3 - ^ ~ Q + r f! = !-18i8Si8^ (-Pt : 2.14 ± 0.63). (4.1) 
B(B — > aJK ) 

In above ratios the hadronic uncertainties are mainly governed by weak annihilation and spectator 
scattering in Ri, R2 and largely canceled out in R3. It is evident that while the predicted R\ is 
barely consistent with the data within errors, theory does not agree with experiment for R2 and 
i?3. This should be clarified by the improved measurements of these modes in the future. 
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TABLE VI: Branching ratios (in units of 10 -6 ) for the decays B ai(1260)vr, 01 (1260) if, 
&i(1235)-7r and 6i(1235)iT. The theoretical errors correspond to the uncertainties due to variation of 
(i) Gegenbauer moments, decay constants, (ii) quark masses, form factors, and (hi) Xr,pa,h, <Pa,h, 
respectively. Other model predictions are also presented here for comparison. In [14], predictions 
are obtained for two different sets of form factors, denoted by I and II, respectively, corresponding 
to the mixing angles 9k 1 = 32° and 58° (see the text for more details). 



Mode 




CMV [15] 


LNP(I)[14J 


LNP(II) 


This work 


Expt. [2 , 3, 4, 5, 6, 7, i0] 






74.3 


4.7 


11.8 


Q 1+0.2+2.2+1.7 
J - ± - 0.2-1. 8-1.1 


12.2 ±4.5 a 


B°^ 


a^7r + 


36.7 


11.1 


12.3 


9 o ,(+2.3+6.2+1.9 
^•3-^-2. 2-5.5-1. 3 


21.0 ±5.4 a 


B°^ 




111.0 


15.8 


24.1 


oo cr+2.5+8.4+3.6 
oz -" J -2.4-7.3-2.4 


31.7 ±3.7 6 


B -» 


a?vr- 


43.2 


3.9 


8.8 


7 fi +0.3+1.7+1.4 
' -°-0.3-l. 3-1.0 


20.4 ±4.7 ±3.4 


B°^ 


a?vr° 


0.27 


1.1 


1.7 


n Q+0.1+0.3+0.7 
u - y -0.1-0.2-0.3 




B -» 




13.6 


4.8 


10.6 


14 4+1.4+3.5+2.1 
^•^-1.3-3.2-1.9 


26.4 ±5.4 ±4.1 




afif" 


72.2 


1.6 


4.1 


IS q+1. 0+14.2+21.1 
lo -°-1.0- 7.2- 7.5 


16.3 ±2.9 ±2.3 


B°^ 


a?K° 


42.3 


0.5 


2.5 


a q+0.3+6.1+9.5 
D - y -0.3-2. 9-3.2 




B- -» 




84.1 


2.0 


5.2 


91 p+\. 2+16.5+23. 6 
zl -°-l.l- 8.5-11.9 


34.9 ±5.0 ±4.4 


i ) — 

B — » 


7^— 


A O A 

43.4 


1 A 

1.4 


2.8 


1 3 Q+0.9+9.5+12.9 
ic5 - y -0.9-5.1- 4.9 




B~^ — > 


b+lT- 


36.2 


6.9 


0.7 


1 i 9+0.3+2.8+2.2 
li - z -0.3-2.4-1.9 




5°^ 


677T+ 


4.4 


« 


« 


n q+0.1+0. 1+0.3 

u - a -0.0-0. 1-0.1 








40.6 


6.9 


0.7 


H A +0.4+2.9+2.5 
ii -^-0.3-2.5-2.0 


10.9 ± 1.2 ±0.9 






0.15 


0.5 


0.01 


i 1+0.2+0.1+0.2 
L - L -0.2-0. 1-0.2 




5- - 




4.2 


4.8 


0.5 


n ^+0.0+0.2+0.4 
u -^-0.0-0. 1-0.2 




B- -» 


b^iT- 


18.6 


4.5 


0.4 


n r+0.3+1. 6+2.5 
J -°-0.3-1.6-1.5 


6.7 ± 1.7 ± 1.0 


B° — >■ 


b\K- 


35.7 


2.4 


0.2 


19 1+1.0+9.7+12.3 
iz - i -0.9-4.9-30.2 


7.4 ± 1.0 ± 1.0 


23° — > 


b°jf 


19.3 


4.1 


0.4 


7 q+0.5+5.4+6.7 
' -°-0.5-2.8-6.5 




£?~ -» 


b-jf 


41.5 


3.0 


0.3 


1 4 n +1.3+ll. 5+13.9 
1^.U_ 1 2 _ 5 9 _ g 3 




B -» 


b\K- 


18.1 


2.6 


0.07 


a 9 +0.5+5. 0+6.4 
D - z -0.5-2.5-5.2 


9.1 ± 1.7 ± 1.0 



a BaBar data only [4[ . 

6 thc average of BaBar [2| and Belle [3| data. 



While the tree-dominated aitr modes have similar rates as pn ones, the penguin-dominated a\K 
modes resemble much more to nK than pK , as first pointed out in [16]. One can see from Eqs. 
(|A3p and ()A7p that the dominant penguin coefficients oO±(a\K) and oFJ^itK) are constructive in 04 
and ag penguin coefficients: 

al( ai K) = a p 4 ( ai K) + r^( ai K), o^ttA") = a^K) + T*<%(irK), (4.2) 

whereas o^pif) = a\{pK)—r^ a p ^{pK) 39[]. Consequently, when the weak annihilation contribution 
is small, B — > a\K and B — > ttK decays should have similar rates. However, if weak annihilation 
is important, then it will contribute more to the a\K mode than the ttK one due to the fact 
that f ai S> f n , recalling that the weak annihilation amplitude is proportional to /b/mi/m 2 &i- By 
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comparing Table IVII of the present work with Table 2 of 39( , we see that the default results for 
the branching ratios of B — > a\K and B —* ttK decays are indeed similar, while the hadronic 
uncertainties arising from weak annihilation are bigger in the former. 

2. B biir,b\K decays 

As for B — ► &i(1235)7r decays, there is a good agreement between theory and experiment. Notice 
that it is naively expected that the mode is highly suppressed relative to the bfir~ one as 

the decay amplitude of the former has the form aiF^ 77 * /fti^Jj 1 (°i being the effective Wilson 
coefficient for the color-allowed tree amplitude) and the decay constant fb ± vanishes in the isospin 
limit. As noted in passing, the LCDA $*Mu) given by (|2.67p is finite even if = 0. This is because 

the coefficient fi^ = \/a\ bl in the wave function of b\ will become divergent if fb x = 0, but the 
combination fb^bx is finite. More precisely, fb^bi is equal to the transverse decay constant of 
the b\ [cf. Eq. (|2.65p ]. Therefore, ^jj 1 (it) can be recast to the form of Eq. (|2.63p which amounts 
to replacing by fj^ in the calculation. Now, one may wonder how to see the suppression of 
6^7r + relative to b\ir~ ? The key point is the term J $¥(x)dx appearing in the expression for the 
effective parameter at [see Eq. (|3.1|) ]. This term vanishes for the b\ meson in the isospin limit. As 
a result, the parameter a\ for the decay 5° — ► 6j~7r + vanishes in the absence of vertex, penguin 
and spectator corrections. On the contrary, a% = c\ + y + • • ■ for the channel B° — > bfir~. This 
explains the suppression of B° — > b^ir + relative to 6^7r~ . After all, the 6j~7r + mode does not 
evade the decay constant suppression. It does receive contributions from vertex and hard spectator 
corrections and weak annihilation, but they are all suppressed. The BaBar measurement of charge- 
flavor asymmetry AC implies the ratio T(B° b{-K + )/T(B° bfir^) = -0.01 ± 0.12 [5]. This 
confirms the expected suppression. 

Since B — > b\K decays receive sizable annihilation contributions, their rates are sensitive to the 
interference between penguin and annihilation terms. As a consequence, the measured branching 
ratios of B — > b\K would provide useful information on the sign of the B — > &i(1235) transition 
form factors. We found that if the form factor V 1 is of the same sign as V<f a \ B(B -► hK) will 
be enhanced by a factor of 2 ~ 3, for example, B(B° — > bfK~) ~ 21 x 10~ 6 which is too large 



compared to the experimental value of (7.4 ± 1.4) X 10 6 [5|. This means that the interference 
between penguin and annihilation contributions should be destructive and the form factors V^ bl 



and Vq 8 " 1 must be of opposite signs. 



We also found that the naive relation B(B~ -» b^K') ~ -» b^R-) holds in QCDF. 

More precisely, QCDF predicts 

R * 55 1% ^ulTl = OM -oT-oT-om (expt : 1.23 ± 0.36), (4.3) 
B(B — > b[K ) 

where the hadronic uncertainty in arises almost entirely from weak annihilation (contribution 
from spectator scattering is negligible). This indicates that the data of b\K~ and b\K~ can be 
simultaneously explained only if the weak annihilation mechanism plays a dominant role in these 
decays. 
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TABLE VII: Same as Table EI except for the decays B 1?i(1270)tt, K x (1270) K, Ki(1400)vr 
and l?i (1400)1? for two different mixing angles 6k 1 = —37° and —58° (in parentheses). In the 
framework of [IJ], only the l?|~(1400)7r° and /C 1 (1400)-7r° modes depend on the mixing angle 9k 1 - 
Note that the results of [LJ] and [15J] shown in the table are obtained for 6k 1 = 32° and 58° (in 
parentheses). 



Mode 



[15] 



M 



This work 



Expt. [8] 



w 

Tf 

B~ 
B~ 
5° 
5° 
B~ 
B~ 



w 
if 

B~ 
B~ 
if 
Tf 
if 
if 
B~ 
B~ 
B° 
if 



Kf(1270)7r+ 
Z?(1270)vr° 
>Z?(1270)7r- 

> ifr(1270)7T° 
ETf (1400)vr+ 
K^(1400)7T° 

>^(1400)vr- 

> iff (1400)tt° 



4.3 (4.3) 
2.3 (2.1) 
4.7 (4.7) 
2.5 (1.6) 
2.3 (2.3) 
1.7 (1.6) 
2.5 (2.5) 
0.7 (0.6) 



7.6 
0.4 
5.8 
4.9 
4.0 



3.0 (1.7) 2 
3.0 6 
1.0 (1.4) 3 



n+0.8+1 

u -0.6-0 

n+o.o+o 

u -0.0-0 
c+0.1+1 
°-0.1-l 
7+0.1+1 
' -0.1-0 
+1.1+1 
1.0-1 
+0.3+0 
0.3-0 
+1.0+2 
0.9-1 



+0.4+1 
-0.4-0 



5+4.2 
9-1.4 
6+1.7 
3-0.6 
8+5.1 
1-1.9 
1+3.1 
7-1.0 
7+9.9 
3-2.8 
7+5.5 
6-1.7 
0+11.6 

6- 3.6 
1+5.2 

7- 1.3 



(3.0j 



/ 7 +0.6+1.3+4.4^ 
\ A - '-0.5-0.8-1.57 
/ n o+O.l+0. 5+1. 7\ 
V u -°-0.1-0. 3-0.6/ 
+0.2+0. 1+2. 7s 
0.2-0. 2-2.2J 
(O c+0.1+1.0+3.2s 
l z -°-0.1-0. 7-1.0/ 
r 9+1. 1+0. 7+2. 6\ 
l z - z -0.8-0.6-1.3/ 
/1 r+0.4+0.3+1.7s 
\ L - o -0.3-0.3-0.9J 
fO O+1.0+0.9+3.0 \ 
l z -°-0. 8-0.9-1. 71 
+0.4+0. 4+1. 2s 
-0.3-0.4-0.5 J 



12.0 ±3.1 



+9.3 
-4.5 



< 25.2 



16.7 ±2.6 



+3.5 
-5.0 



< 21.8 



(1.0 



l?f (1270)1?+ 
1?+(1270)1?- 

> 1?°(1270)1?- 

> Iff (1270)1?° 
1? (1270)1? 
1?°(1270)1? 
l?f (1400)1?+ 
1?+ (1400)1?- 

> l?°(1400)i?- 

> l?f (1400)1?° 
1?°(1400)1? 
1?°( 1400)7?° 



0.22 (0.22) 
0.02 (0.75) 
0.02 (0.70) 
0.20 (0.20) 



0.12 (0.12) 
4.4 (3.9) 
4.1 (3.6) 
0.11 (0.11) 



ho. 
-0. 
+0. 
-0. 
+0. 
-0. 
+0. 

0. 
+0. 
0. 

241°; 



01 
06 
25 
05 
30 



011° 

081° 
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-0 



01+0. 

00- 0. 
01+0. 

01- 0. 
01+0. 

01- 0. 
02+0. 

02- 0. 

16+1. 
15-0. 
01+0. 
01-0. 
01+0. 
01-0. 
00+0. 
00-0. 
08+0. 
08-0. 
00+0. 

00- 0. 
01+0. 

01- 0. 
08+0. 
07-0. 



00+0.02 
00-0.01 
00+0.46 
00-0.06 
15+0.39 
08-0.09 
07+0.10 
03-0.04 
13+1.43 
61-0.61 
11+0.33 
07-0.13 
00+0.23 
00-0.09 
00+0.04 
00-0.00 
15+0.81 
12-0.26 
01+0.14 
00-0.01 
17+0.59 
06-0.08 
13+0.92 
11-0.32 



0.01 



hO.OO 
-0.00 

n n4+ - 01 

U - U4 -0.01 

n 22+°- 01 

u - zz -0.01 



0.05 
2.10 



0.02 
0.01 

+0.13 
0.13 



26+ - 10 
u - ZD -o.oi 

07+ - 02 

u - u '-0.02 

n ni+ - 00 
u - Ui -o.oo 



0.22 



-0.07 
-0.07 

01+°+° 

U.U1_ _q 

n in+ - 02 

u - iu -0.02 
u - zo -0.07 



0.00 
-0.00 
+0.00 
-0.00 
+0.12 
-0.07 
0.09 
0.03 
1.23 
0.65 
0.12 
0.08 
0.00 
0.00 
0.00 
-0.00 
0.07 
-0.07 
02+0 
00-0 
fO.21 
-0.08 
+0.08 
-0.07 



0.02 s 

-0.01/ 
0.27 s 
0.04 J 
0.39 \ 
0.12/ 
0.10s 
0.04/ 
1.31s 
0.57^ 
0.47 s 
0.17/ 
0.16s 
0.06/ 
0.16s 

-0.00/ 
0.24 s 

-0.13/ 

04 s 

00 ) 

f 0.15 s 
-0.10 J 
f 0.31s 
-0.15/ 



3. B ->l?i (1270) (tt, 1?),1?i (1400) (tt,1?) decays 



It is evident from Table IVlTl that the central values of the calculated branching ratios in QCDF 
for l?j~(1270)7r + and l?-f (1400)7r + are too small compared to experiment. This is not surprising 
as the same phenomenon also occurs in the penguin dominated B — ► PV and B — ► VV decays. 
For example, the default results for the branching fractions of B — ► K*n obtained in QCDF are 
in general too small by a factor of 2 ~ 3 compared to the data [39(] . This suggests the importance 
of power corrections due to the non-vanishing pa and ph parameters or due to possible final-state 
rescattering effects from charm intermediate states 42|. It has been demonstrated in 39J that in 
the so-called "S4" scenario with pa = 1 and non- vanishing 6 a , the global results for the VP modes 
agree better with the data. It has also been shown in 43] that the choice of pa^ a = 0.6e _l4 °° 
will allow one to explain the polarization effects observed in various B — > VV decays. While large 
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TABLE VIII: Same as Table IVTl except for the decays B — > fin, f\K, h\ir and h\K with fi = 
/i(1285), /i(1420) and h\ = /ii(1170), _i(1380). We use two different sets of mixing angles, namely, 
3Pl = 27.9° and 9i Pl = 25.2° (top), corresponding to 9 Kl = -37°, and 3Pl = 53.2°, 9i Pl = -18.1° 
(bottom), corresponding to 9k 1 = —58°. 



Mode 



Theory 



Mode 



Theory 



B~ 



> /i(1285)tt- 
/i(1285)7r° 

5- -> /i(1285)K- 
/i(1285)__° 

> /ii(1170)vr- 

%(1170)7T° 

5- -> hi(U70)K- 
B° -> /ii(1170)__° 



5° 
£T 
5° 



0.26 



r n+0.3+1.3+0.7 
d - z -0.2-l. 0-0.2 
+0.03+0.14+0.29 



14i 



-0.03-0.07-0.08 
h3.0+7.5+12.4 
'-2.6-3.9- 5.2 
1zL fi +2. 7+7.5+11. 9 
ly: - o -2.3-3.9-5.0 
/i n+0.4+0. 9+0.8 
^■°-0.3-0. 8-0.7 



B- 
B { 



0.19 



+0.06+0.05+0.07 
-0.04-0.03-0.01 



10.1 
10.1 



+4.7+2.1+7.3 



-3.1-1.4- 



.1 



+4.2+2.2+7.2 



-2.8-1.5- 



.1 



> /i(1420)tt- 
/i(1420)tt 

B~ -> /i(1420)A- 
/i(1420)__° 

> /ii(1380)tt- 

/ii(1380)7r° 

> /ii(1380)A- 
U° -> /ii(1380)K° 



n nK+O-Oi+O-Oi+O- 00 
u - uu -o.oo-o.oo-o.oo 



0.003; 



6.0] 



0.006 



-0.003+0.002+0.003 
-0.002-0.001-0.002 
+1.7+1.9+9.0 
1.5-1.3-3.1 
r r+1. 6+1.8+8.4 
J - J -1.3-1.2-2.8 
n i 7+0.03+0.06+0.04 
U- 1 ' -0.02-0.06-0.03 
+0.009+0.005+0.007 
-0.004-0.004-0.002 
i 9 7+7.1+9.2+211.4 
-5.1-4.7- 10.8 
i i q+6. 4+8.5+188. 5 
ll--J_4.6-4.3- 9.6 



B~ 



► /i(1285)tt- 
/i(1285)7t° 

B /i(1285)i_" 
3° -> /i(1285)i_° 

► /ii(1170)vr- 

/ll(1170)7T 

f?- /ii(1170)i_- 

if -» /ii(ii7o)_r° 



. fi +0.2+l. 1+0.6 
^■"-0.2-0.9-0.2 



0.20 



+0.02+0.12+0.24 
-0.02-0.06-0.06 



5.2 



+0.9+3.1+ 9.1 



0.8-1.5-10.0 

5 2+0.8+3.2+3.0 



- -0.7-1.5-1.4 
2+0.3+0.3+0.3 
'-0.2-0.3-0.3 
n 1fi +0.08+0.0l+0.06 
u - ±o -0.05-0.01-0.04 
-5.8+1.9+23.0 
-3.7-1.1- 8.2 
+5.3+1.9+21.1 
-3.4-1.2- 7.7 



1, 



11.3] 
10.9 



> /i(1420)vr- 
/i(1420)tt° 

> /i(1420)i_" 
/i(1420)__° 

> fci(1380)ir - 
/ii(1380)7r° 

_? _ -► /»i(1380)__ _ 
1?° -» /ii(1380)__° 



B~ 
B° 
B- 
B° 
Br 



n c:n+0.06+0 
u - jy -0.05-0 



0.05 



+0.02+0 
-0.01-0 



13.8 



+4.0- 
-3.3- 



iq 1+3.7 + 

+0.2 
-0.1 



2.9 



0.04; 



-0.00+0 
-0.00-0 

5 6 +0 - 9 " 

O.D_q 7 _ 

-+0.7- 
-0.7- 



5.5] 



.18+0.10 
.13-0.05 
03+0.04 
02-0.02 
5.6+17.1 
3.2- 6.3 
5.4+16.2 
3.1- 5.9 
fO.6+0.4 
-0.6-0.4 
03+0.04 
02-0.01 
-2.3+1.4 
-1.2-1.9 
-2.4+1.5 
-1.2-2.0 



power corrections from weak annihilation seem to be inevitable for explaining the Km rates, one 
issue is that large weak annihilation may destroy the existing good agreement for afK~ and b\K~ 
modes. 

We notice that while i_i(1270)7r rates are insensitive to the mixing angle 9k\i the branching 
fractions of i_i(1400)7r are smaller for 9k 1 = —58° than that for 9k\ = —37° by a factor of 2 ~ 3 
due to the dependence of the i_i(1400) decay constant on 9k 1 , recalling that /ki(i4oo) ~ H2 (235) 
MeV for 9 Kl = -58° (-37°) [cf. Eq. £T32l) ]. The current measurement of B° -» i_f(1400)vr + 
favors a mixing angle of —37° over —58°. 

Just as the case of B —* KK* decays, we find that the branching ratios of B — * i_i(1270)i_ 
and i_i(1400)i_ modes are of order 10 -7 — 10 -8 except for the decay B° — » k\{127Q)K° which can 
have a branching ratio of order 2.3 x 10~ 6 . The decay modes K±K + and K^K~ are of particular 
interest as they are the only AP modes which receive contributions solely from weak annihilation. 



4- B — ► /i(7r, K), hi(Tr, K) decays 



Branching ratios for the decays B — > fin, f\K, h\n and h\K with f\ = f\ (1285), f\ (1420) and 
hi = hi (1170), h\ (1380) are shown in Table IVTril for two different sets of mixing angles: (i) 9a P = 
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0° 90° 180° 

FIG. 1: Branching ratios of B~ -> /i(1285)if- (solid line) and B~ -> f 1 {U20)K~ (dashed line) 
versus the mixing angle 63 p, . The physical mixing angle is either 28° or 53°. For simplicity, only the 
central values are shown here. The horizontal line is the experimental limit on B~ — > fi(1285)K~ . 

53.2° and 9i Pl = -18.1°, corresponding to 9 Kl = -37°, and (ii) 9 3Pl = 27.9° and 9i Pl = 25.2°, 
corresponding to 9k 1 = —58° (see Sec. II. A). 6 Their branching ratios are naively expected to be 
of order 10 -6 ~ 10 -5 except for the color-suppressed fin and h\n° modes which are suppressed 
relative to the color-allowed one such as /i(1285)-7r~ by a factor of |«2 / Q-i | 2 /2 ~ 0(0.03 — 0.08). 
However, an inspection of Table I Villi shows some exceptions, for example, B(B~ —* t /i(1420)-7r~) <C 
B(B- -> /i(1285)vr-) for both sets of 9 3Pl and B(B~ -> /ii (1380)0 < B{B' -> fri (1170)0 for 
9i Pl = 25.2°. These can be understood as a consequence of interference. The decay amplitudes for 
the tree-dominated channels h\n~ are given by 

A(B~ -> fci (1380)0 oc -V Bhl sm8i Pl + V Bha cos0i Pl , 

A(B~ -> fci (1170)0 oc V^ 1 cos0i Pl + V^ 8 sin0i Pl . (4.4) 

Since the form factors V Bhl and V Bhs are of the same signs (cf. Table ITTT|) . it is clear that the 
interference is constructive (destructive) in the hi(1170)n~ mode, but destructive (constructive) 
in fti (1380)71- for 9i Pl = 25.2° (-18.1°). This explains why B(B~ -> /n (1380)0 < B(-B~ -» 
fa (1170)0 for 6»i Pl = 25.2° and B(B" -» /*i(1380)O > £(£~ -> fa (1170)0 for 6>i Fl = -18.1°. 
Therefore, a measurement of the ratio R§ = B(B~ — > /ii(1380)7i— )/jS(-B~ — > /ii(1170)7r~) will help 
determine the mixing angle 9ip r Likewise, information on the angle 93 Pl can be inferred from the 
ratio R 6 = B(B -> fx{lA3Q)K)/B(E -> /i(1285)K): # 6 > 1 for # 3Pl = 53.2° and #6 < 1 for 
0s Pl = 27.9°. 

The preliminary BaBar results are 0] 

B(B- -> /i(1285)iT-)B(/i(1285) -» tjtttt) < 0.8 x 10~ 6 , 
-> /i(1420)iT~)fi(/i(1420) -► t/tttt) < 2.9 x 10~ 6 , 
B{B~ -> /i(142O)if-)0(/i(142O) -► ICgA^O < 4.1 x 10" 6 . (4.5) 



6 There are predictions for the decay rates of B — > fin, f\K, hin and /iiif in [15|. Since the /1 and hi 
states are not specified there, we will not include them in Table fVIIII for comparison. 
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TABLE IX: Direct CP asymmetries (in %) in the decays B ai(1260)vr, ai(1260)K, &i(1235)vr 
and b\{122>b)K. See Table IVI1 for the explanation of theoretical errors. Experiments results are 
taken from 



Mode 



Theory 



Expt. 



W 

if 

5° 

B- 
B- 



Mode 



Theory 



Expt. 



w 

B° 
B~ 
B~ 



w 
if 

B- 
B- 



a{ 7T 
af 7r + 
a?vr° 

>a?7r- 



-3.6 
-1.9 
60.1 
0.5 
-4.3 



0.1+0.3+20.8 
0.1-0.5-20.2 

+0.0+0.0+14.6 
0.0-0.0-14.3 

+4.6+6.8+37.6 
4.9-8.3-60.7 
0.3+0.6+12.0 
0.2-0.3-11.0 

+0.3+1.4+14.1 
0.3-2.2-14.5 



7 ±21 ±15 

15 ± 15 ± 7 



* a 1 K 
► a\K- 



2.6 
-7.7 



n 8+°-° 



0.0+0.7+10.1 
0.1-0.7-11.0 
+0.6+2.1+6.8 
0.6-2.2-7.0 
0.1+0.6 
0.1-0.0 



-16 ± 12 ± 1 
12 ± 11 ±2 



o A +0.3+1.4+10.3 
"■^-0.3-1.6-12.0 



b{7T- 
b\^ 



-4.0 
66.1 



+0.2+0.4+26.2 
0.0-0.6-25.5 

+1.2+7.4+30.3 
1.4-4.8-96.6 



co 4 +6.4+9.0+5.2 
"-> -^-6. 3-7.3-4.7 



-36.5 



+4.4H 
-4.3- 



■18.4- 
17.7- 



-82.2 
59.6 



B 
B° 
B~ 
B~ 



bJK~ 

bW° 
>b^K° 
>b\K~ 



5.5 
-8.6 
1.4 
18.7 



0.2+1.2+47.2 
0.3-1.2-30.2 
0.8+3.3+ 8.3 
0.8-4.2-25.4 
+0.1+0.1+5.6 
0.1-0.1-0.1 
+1.6+7.8+57.7 
1.7-6.1-44.9 



-7 ±12 ±2 



-46 ± 20 ± 2 



0.9 



+0.6+2.3+18.0 
-0.4-2.7-20.5 



5 ± 16 ± 2 



a taken from 



44|. 



Since B(/i(1285) -> 7?7rvr) = 0.52 ±0.16 fl, the upper limit on B{B- /i(1285)if-) is inferred to 
be of order 2.0 x 10 -6 . However, we cannot extract the upper bound for the /i(1420)iT _ mode due 
to the lack of information on £(/i(1420) -> t/tttt) and 23(/i(1420) -> Ks^if?). In Fig. Q]we plot 
the branching ratios of B~ —* /i(1285)A'~ and B~ — > /i(1420)i"C~ as a function of Ozp x . We see 
that the branching fraction of the former at 9-z p 1 = 53° is barely consistent with the experimental 
limit when the theoretical errors are taken into account. Note that the mixing angle dependence 
of the /i(1420)i"C~ mode is opposite to that of f\(1285)K~ . At this moment, it is too early to 
draw any conclusions from the data. Certainly, we have to await more measurements to test our 
predictions. 



B. CP asymmetries 

Direct CP asymmetries for various B — > AP decays are summarized in Tables IIXIIXII Due to the 
large suppression of B° — > 6^7r + relative to B° — > b^ir~ , direct CP violation in the later should be 
close to the charge asymmetry Ab l7T defined below in Eq. (|4.1ip which has been measured by BaBar 
to be —0.05 ± 0.10 ± 0.02 [5J. The default results for direct CP violation vanishes in the decays 
B° — > KfK~ and B° — > K + (see Table |X|) as they proceed only through weak annihilation. 
The major uncertainty with direct CP violation comes from the strong phases which are needed 
to induce partial rate CP asymmetries. For penguin dominated decays, one of the main sources of 
strong phases comes from defined in Eq. (|3. 19j) which is originated from soft gluon interactions. 
It is nonperturbative in nature and hence not calculable. 

The experimental determination of direct CP asymmetries for afir^ and a^7r + is more com- 
plicated as B° — > afir^ is not a CP eigenstate. The time-dependent CP asymmetries are given 



27 



TABLE X: Direct CP asymmetries (in %) the decays B -► i^i(1270) vr, K 1 (1270)K, Ki(1400)vr 



and Ki(lAOO)K for two different mixing angles 9k 1 = —37° (top) and —58° (bottom). 



Mode 


Theory 


Mode 


Theory 


B° 


-► iff (1270)vr+ 


oo c+2. 8+10. 0+26.0 
oo - o -3.3-13.0-42.9 


B° 


X 1 + (1400)tt+ 


14 q+9.0+1. 7+45. 3 
1 ^-°-9.5-1.0-48.8 


B° 


^ir?(1270)vr° 


q 9 r+1.6+7.4+18.7 
oz -°-1.8-7.6-14.8 


B° 


->■ K°(1400)vr 


1 r+0. 2+0.9+2.3 
J-- -!. 2-0.8-2. 2 


B~ 


^^(1270)vr- 


n O+0.3+0.1+3.3 
—u -°-0. 3-0.2-4.2 


B~ 


->■ K°(1400)vr- 


9 9+0.4+0.1+1.5 
z - z -0. 3-0.1-1. 2 


B° 


-► iff (1270)tt° 


oo q+1.9+7.1+24.5 
J0 -°-l. 8-9.1- 2.5 


B° 


(1400)vr 


,9 q+7.3+2.2+40.4 
— ±z.y_ 6 5 _ 2 7 _ 9 x 


5° 


-► Kf(1270)if+ 


n+0+0+7.3 
u -0-0-7.3 


B° 


-► iff (1400)^+ 


n +0+0+13.5 
u -0-0-13.5 


5° 


-► K+(1270)if- 


H+0+0+40.9 
u -0-0-40.9 


B° 


-► iq^Moo)^- 


n +0+0+88.9 
u -0-0-88.9 


B° 


^(1270)^° 


Qi n +3.5+4.4+10.8 
° J - u -3.1-3.5-11.2 


B° 


-► Z?(1400)K° 


7 1 +0.7+14.9+51.1 
'■ L -0.7- 5.1-14.0 


5° 


-► #?(1270)Z° 


17 q+6.0+0.9+3.9 
1 ' - y -6.4-0.8-3.0 


B° 


-► ^(MOO)^ 


ct; c+4.0+1.4+12.9 
od - d -4.0-l. 4-18.0 


B~ 


-»• Kl(127Q)K- 


17 n+5.8+4.2+65.0 
i '- z -5.4-2.9- 7.6 


B~ 


-► if?(1420)if- 


n+7.7+1.5+23.7 
^°- u -7.5-1.0-10.8 


B~ 


-► iff(1270)if° 


A n 9+6.7+2.1+131.3 
^ u - z -8. 3-5.0- 11.7 


B~ 


-► K 1 ~(1420)i ; s: 


1 7 q+22.1+26.4+101.8 
i '-°-22.1-36.2- 31.2 


B° 


-► iff (1270)vr+ 


qq R +2.6+ 8.5+31.2 
00 -°-2. 3-10. 1-50.7 


B° 


iff (1400)vr+ 


qq 9+8.O+1. 5+40. 7 
oy - z -2.9-2.9-35.4 


B° 


^Zi(1270)vr° 


9Q r+1.4+6.8+19.7 
zy - D -1.4-7.9-23.5 


B° 


->■ K°(1400)vr 


n 1+4.3+1.3+6.0 
u - ± -5.0-1.2-5.6 


B~ 


^^(1270)^- 


n c+0. 2+0.0+2.7 
u -°-0. 2-0.2-2. 2 


B~ 


->■ i?°(1400)vr- 


q 1 +0.2+0.9+3.0 
o - 1 -0.1-0.9-1.7 


B° 


-► (1270)vr° 


nn „j-n 5+5 1+26 5 
JZ -°-0.5-6.7- 0.7 


B° 


-»• Xi"(1400)7T° 


^r. n+7 2+5 6+33 9 
— ^ z - z -8.5-8. 7-12.1 


B° 


-► iff (1270)K+ 


n +0+0+27.0 
u -0-0-27.0 


B° 


-► Xi"(1400)K + 


n +0+0+19.7 
u -0-0-19.7 


B° 


-► K?(1270)K- 


n+o+o+40.1 

u -0-0-40.1 


B° 


-► Jff (1400).K' _ 


n +0+0+48.1 
u -0-0-48.1 


B° 


^Zi(1270)K° 


is n+°- 5 + 1 - 2 + 1 - 3 

io - u -0. 6-0.7-0.1 


B° 


-► Zi(1400)K° 


0A Q+2.2+2.4+10.1 
z ^-°-2.0-5.2-56.3 


B° 


-► if 1 °(1270)K° 


n n +i. 0+2. 4+13.0 

— 1.2—2.1—11.0 


B° 


-► Kf(U00)K° 


R q n+4.3+1. 7+19.3 
oo - u -3.0-2.4-23.8 


B- 


-»• K?(1270)K- 


O .+3.1+3.5+41.7 
y -^-3. 3-2.5- 3.3 


B- 


-► ^(1420)^ 


^Q. c+2.4+1.2+43.3 
oy -°-2.8-1.4- 3.9 


B- 


-► iff (1270)K° 


qq ^+6.7+4.3+127.7 


B~ 


-► iff(1420)if° 


R i 1+21.6+23.2+36.3 
u± - 1 -18.0-41.8-25.5 



by 

r(B°(t) of vr^) - T{B°{t) of ^ 



r(B°(i) -f of TFT) + r(B°(t) -> of 7T+) 

= (S± AS)sm(Amt) - (C± AC)cos(Amt), (4.6) 

where Am is the mass difference of the two neutral B eigenstates, S 1 is referred to as mixing-induced 
CP asymmetry and C is the direct CP asymmetry, while AS and AC are CP-conserving quantities. 
Defining 

A + _ = A(B° a+7T) , A_ + = ^(B° orvr + ) , 

= A(B° o7 7T+) , i+_ = A(B° of O, ( 4 - 7 ) 

and 

A - Q A+ - \ - q A ~ + (A 

where q/p = e~ 2t @ for a\7r modes, we have 

C l AC _ l-|A+-| 2 _ l^+-| 2 -|i+-| 2 c Ac _ l-|A-+I 2 _ |^-+| 2 -|A- + | 2 
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TABLE XI: Direct CP asymmetries (in %) in the decays B — ► /i7r, /i-ftT, fai7r and foil? with 
/j = f x (1285), /i (1420) and h x = fa (1170), fa (1380). We use two different sets of mixing angles, 
namely, 9sp 1 = 27.9° and 9ip 1 = 25.2° (top), corresponding to 9k 1 = —37°, and 03p 1 = 53.2° and 
8i Pl = —18.1° (bottom), corresponding to 6k 1 = —58°. 



Mode 



Theory 



Mode 



Theory 



B~ 
B° 



► /i(1285)vr- 
/l(1285)7T° 

B~ -» /i(1285)iT- 
/i(1285)Z° 

► fcx(1170)7r~ 

/»l(1170)7T° 

B -» /ii(1170)if- 
B° -> fa (1170)7?° 



5° 
J3- 
B° 



7 q+0.4+0.5+28.0 
' - J -0.5-0.6-27.5 
1 3 0+1-7+3.2+52.6 
10 -°-l. 7-3.8-58.7 
9 c:+0.3+0.8+6.7 
z -°-0.2-0.7-8.0 
0.1+0.4+2.0 
0.1-0.4-2.6 
+1.1+1.8+3.4 
-1.1-2.3-3.7 
7.5+6.5+58.8 
-6.8-7.5-76.7 
o c+0.8+0.3+14.1 
°-°-0. 7-0.2-16.5 
i y+0.2+0.2+2.9 
L - ' -0.1-0.1-2.7 



B~ 
B° 



1.9 
-11.1 
31.6 



> /i(1420)7r- 
/i(1420)7r° 

B- -> A (1420)1?" 
/i (1420)1?° 

► /ii(1380)vr- 
fa(1380)vr 

B- -» fa(1380)l?- 



B° 
B° 



B° 



fa(i380)i?° 



-4.1 



+0.9+0.6+44.6 



-1.1-0.7-44.2 
n A n+10.3+8.0+112.2 
-C<±.U_ g.i_3.3_ 66.4 

n 0+0.1+0.1+1.6 

u -°-0. 1-0.1-1. 7 



1.0 



+0.1+0.2+0.9 
-0.2-0.2-0.8 
i o 9+2.1+3. 6+23.3 
i0 - z -2.3-5.4-23.3 
no 7+13.9+13.9+124.1 
— 18.7 — 9.7- 72.4 



1.3 



+0.2+0.8+12.; 



0.5-0.6-18.2 
1 O+0.3+0.4+0.4 
i -°-0.3-0.5-3.7 



B" 
B° 



► /i(1285)vr- 
/i(1285)tt° 

B~ -» /i(1285)K- 
A(1285)F° 

► fa(1170)7T 
fa(1170)7r° 

B~ /n (1170)1?- 



B° 
B~ 
B° 



B° 



fa(1170)l?° 



-7.1 
14.7' 



+0.4+0.5+28.6 
-0.4-0.6-28.0 
-1.7+3.0+57.0 



-1.9-3.8-64.0 
-0.3+0.9+3.1 
-0.3-0.7-1.4 
9 1+0.2+0.6+2.5 
jS - i --0.1--0.5-3.1 
o .+0.7+1.3+7.2 
~ o ' 4 -0.8-1.7-7.9 
-7.2+2.8+46.5 
-5.6-3.2-52.1 
-2.2+0.8+13.7 
-1.4-0.9-12.4 
n 7+0.2+0.2+2.1 
u -' -0.2-0.2-2.1 



B- 
B° 



2.5; 



24.9; 
5.0 



> /i(1420)tt- 
/i(1420)7r° 

B- -» A (1420)1?" 
/i(1420)l? 

► /ii(1380)vr- 
fa(1380)7r° 

B- fa (1380)1?" 



B° 
B" 
B° 



B° 



fa (1380)1?° 



-3.8; 

21.0 



-0.3+0.4+26.4 
-0.4-0.4-26.0 
+3.2+7.0+40.2 
2.7-6.5-45.0 
i q+0. 3+0.3+3.4 
1 -°-0. 2-0.3-3. 5 
i 1+0.2+0.2+1.1 
1 - 1 -0. 1-0.2-1. 2 
io q+1. 4+2.3+1.4 
io -°-1.3-2.9-1.4 
ra A+3-7+11.2+ 53.5 
JU - U -4.1-13.6-116.3 

i +0.4+3.6+4.4 
°- L -0. 2-3.5-4.0 

1 E+0.5+0.3+1.6 
1 - d -0. 6-0.5-1. 5 



S-AS ^ ■ u ,2 = WTU .7 - 



and 

2ImA + , 2Im(e" 2 ^A f -,4 + _) 
1 + IA+-I 2 L4 + _| 2 + L4+_| 2 • 
2ImA_ + _ 2Im(e- 2ifS A^ + A*_ 
1 + |A_+| 2 ~ \A_ + \i + \A_ 
Hence we see that AS describes the strong phase difference between the amplitudes contributing 
to B° — > afir^ and AC measures the asymmetry between T(B° — ► a^rr") + T(B° — > a7Tr + ) and 
r(B° -► arvr+) + T(B° -» of tt"). 

Next consider the time- and flavor-integrated charge asymmetry 

^ - l^-l 2 + l^+-l 2 -|-4-+l 2 -|^-+l 2 (411) 



Then, following 38J] one can transform the experimentally motivated CP parameters A ai n and C ai7T 



into the physically motivated choices 



_ \hc-+\ 2 -1 _ \k+-\ 2 -1 

A atn~ = | K -+|2 + 1 > A a-^+ = | K +-|2 + 1' ^ l2 > 



with 

+- q A-+ -+ q A 



(4 - 13> 
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TABLE XII: Various CP parameters for the decays B° — > afir^ (top) and B° — > bfir^ (bottom). 
The parameters S and AS are computed for (3 = 22.1° and 7 = 68.0°. Experimental results are 
taken from M, 15J] . 



Parameter 


Theory 


Experiment 


C 

s 

AC 
AS 


n nn 1 ? +0.001 +0.002 +0.043 
u.uuo n nn9 —0 nrn — n ru^ 

— U.UUi — U.UUO — U-KJ^iO 

n n9 +0.00 +0.00 +0.14 

u - uz -0.00 -0.00 -0.14 
n 07+O.OI +0.05+0.09 

-0.01-0.08-0.16 
n A A +0.03+0.03+0.03 

-0.04 -0.05 -0.04 

n m +0.00+0.00+0.02 

U - Ui -0.00-0.00-0.02 


-0.07 ±0.07 ±0.02 
-0.10 ±0.15 ±0.09 
0.37 ±0.21 ±0.07 
0.26 ±0.15 ±0.07 
-0.14 ±0.21 ±0.06 


1 

«<7fr 

«eff 


l y '- z -0.3-0.6-2.5/ 
ft n7 n +0.5 +3.6 +6.6no 
l iu '- u -0. 5-2.3-3. 7/ 

/ im n+0- 4 + 2 - 3 + 5 -<V 

l iuz - u -0.4-1.5-3.lJ 


(78.6 ± 7.3)° 


c 
s 

AC 
AS 


n n«+0- 01 +0- 01 +0- 23 

U - UD -0.01 -0.01 -0.23 

n n o +0.01 +0.01 +0.06 

u - uo -0.02 -0.02 -0.01 

n n^+0- 03 +o o 2 +0- 15 

u - uo -0.03 -0.02 -0.26 
n qr +0.03 +0.02 +0.08 
u - yu -0.03 -0.03 -0.01 
n -1 9 +0.04+0.04+0.08 
u - iz -0. 03 -0.04 -0.09 


-0.05 ±0.10 ±0.02 
0.22 ±0.23 ±0.05 a 

-1.04 ±0.23 ±0.08 


«rfr 

"eff 


M07 fi+0.7+3.5+155.4\o 
l ±u '-°-0. 2-4.9- 17.8/ 
Hni q+0.4 +2.1 +4.9\o 
V iui -° -0.4-1.4-8.6/ 
nn/l A +0.6+2.6+80.4no 
l iu ^-0.3-2.1- 1.6/ 





°Our definition of C in Eq. ()4.9|) has an opposite sign to that defined in [5[ for B — ► b\ix decays. 



Hence, 

.4 



F(B° -> a+rr-) - T(B° -> cCrC) _ «4 ai7r - C ai?r - AC, 



ai7r 



r(B°^a+7r-)±r( J B ^ar- 7 r+) " 1 - AC ai7 , - Aa 1% C ai7T 
A _ r(^° ^ a^rC) - r(£?° ^ a+rr-) _ ^ + C ai ,r + A ain AC ai7T 
° lf+ r(5°^af7r+)±r( J B°^a+7r-) 1 + AC ai7r + Ai7rC ai7r 

Note that the quantities -A a ±„.T nere correspond to -A a T„.± defined in [38J]. Therefore, direct 
CP asymmetries A a + T _ and -A a -„.+ are determined from the above two equations. 



Defining the effective phases 



1 ^_ 1 



a 



off 



- argK + = - arg (e 2i/3 A- + A*, 



2 2 



"off = 5 ar g K ~ + = \ ar S (e~ 2 ^+-^+) , (4.15) 

which reduce to the unitarity angle a. in the absence of penguin contributions, we have 

l ( + 
"off = 2(«eflf + «cfr 
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S + AS 



arcsm 



+ arcsin 



S- AS 



^/l-(C-ACf 



(4.16) 



This is a measurable quantity which is equal to the weak phase a in the limit of vanishing penguin 
amplitudes. 

Parameters of the time-dependent decay rate asymmetries of B° — > afn^ are shown in Table 
IXIIi It appears that the calculated mixing-induced parameter S is negative and the effective 
unitarity angle a e ff deviates from experiment by around 2a. As pointed out by one of us (K.C.Y.), 
this discrepancy may be resolved by having a larger 7 > 80° (see Fig. 1 of la]). Further precise 
measurements are needed to clarify the discrepancy. For B° — > bfn^ decays, the predicted AC 



T indicates that the B 



■0 



agrees with experiment. The fact that this quantity is very close to 

b^n + mode is highly suppressed relative to the B° — > b^n^ one, recalling that AC here measures 
the asymmetry between T(B° b^n') + T(B° -> 67vr+) and T(B° -> b{n + ) + T(B° 



b+TT-). 



C. Comparison with other works 

There are several papers studying charmless B — > AP decays: Laporta, Nardulli, and Pham 
(LNP) Q (see also Nardulli and Pham [13]), and Calderon, Muhoz and Vera (CMV) Q. Their 
predictions are shown in Tables IVII and IVIII Both are based on naive factorization. While form 
factors are obtained by CMV using the ISGW2 model, LNP use ratios of branching ratios to 
deduce ratios of form factors. Hence, the relevant form factors are determined by factorization and 
experimental data. Specifically, LNP found 



V Bai (0) V Q BKl A (0) 

± Bp (o) * 4! K 

V Bbl {0) V BKlB (0) 



h cos 6k 1 — k sin 0k 1 , (4-17) 



A Bp (0) A™*(0) 

where the kinematic factors h, k are defined in [14]. Therefore, LNP used two different sets of 
form factors corresponding to different mixing angle values Qk\ = 32° and 58°. 7 Since LNP only 
considered factorizable contributions to B — > AP decays, it turns out that in B — > K\ir decays, 
only the 1^(1400)71"° and A^(1400)7r° modes depend on the mixing angle 9k 1 - The other K\T[ 
rates obtained by LNP (see Table lvn|) are mixing angle independent. 

The predicted rates for aiir, bin and b\K modes by CMV are generally too large compared 
to the data, presumably due to too big form factors for B — > a\(bi) transition predicted by the 
ISGW2 model. The relation B(JT -» oJ"7r _ ) > B(W — ► a 1 7r + ) is in conflict with experiment. A 
noticeable result found by CMV is that B(B~ -> K{ (1400)^°) and B(B° K°(1400)i^°) are of 
order 10 -6 , while in QCDF they are highly suppressed, of order 10 -7 — 10~ 8 . 

It is clear from Table IVTl that in the LNP model, the form factors (|4.17j) derived using 6k 1 = 32° 
give a better agreement for a\ii modes, whereas 9k 1 = 58° is preferred by bin and b\K data. This 



7 Since the B — > K\a and B — > K\b form factors obtained in the ISGW2 model are opposite in sign, the 
preferred mixing angle 9k x should be negative, as discussed in Sec. II. A. 
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indicates that the data of cli(tt, K) and b\{ir, K) cannot be simultaneously accounted for by a single 
mixing angle 9k 1 in this model. 

Branching ratios of B — > f±P and B — ► h\P are found to be of order 1CP 5 for P = ir^, 77, 7/ , X 
and 0(1CP 7 ) for i 3 = 7r° by CMV. In general, the CMV's predictions are larger than ours by one 
order of magnitude. 

V. CONCLUSIONS 

In this work we have studied the two-body hadronic decays of B mesons into pseudoscalar 
and axial vector mesons within the framework of QCD factorization. The light-cone distribution 
amplitudes for 3 Pi and 1 P\ axial- vector mesons have been evaluated using the QCD sum rule 
method. Owing to the G-parity, the chiral-even two-parton light-cone distribution amplitudes of 
the 3 Pi ( 1 -Pi) mesons are symmetric (antisymmetric) under the exchange of quark and anti-quark 
momentum fractions in the SU(3) limit. For chiral-odd light-cone distribution amplitudes, it is 
other way around. Our main conclusions are as follows: 

• Using the Gell-Mann-Okubo mass formula and the i^i(1270) and Ki(1400) mixing angle 
9 Kl = -37°(-58°), the mixing angles for 3 P 1 and 1 P 1 states are found to be 9 3Pl ~ 28°(53°) 
for the /i(1420) and /i(1285) and 6i Pl ~ 25° (-18°) for /ii(1170) and the /n(1380), respec- 
tively. 

• The predicted rates for af (1260)^, bf (1235)^, 6?(1235)7r _ , a^K~ and bfK~ modes are 
in good agreement with the data. However, the expected ratios B(B~ — ► aiir~)/B(B° — ► 
afn-) < 1, B(B- -► a{TT°)/B(B° -► a{n+) ~ \ and B(B~ -► b° 1 K-)/B{B° -► b+R-) ~ ± 
are not borne out by experiment. This should be clarified by the improved measurements of 
these decays in the future. 

• One of the salient features of the 1 P± axial vector meson is that its axial- vector decay constant 
is small, vanishing in the SU(3) limit. This feature is confirmed by the observation that 
F(B° — ► b^ir + ) <C T(B° — > 6^7r _ ). By contrast, it is expected that r(P° — > a7/7r + ) S> 
r(P° — > a^7r~) due to the fact that f ai S> /tt- 

• While B ^ a\TT decays have similar rates as that of B — > p7r, the penguin-dominated decays 
P — ► a±K resemble much more to the ttK modes than pK ones. However, the naively 
expected ratio B(B~ -»• a{K°)/B(B° -► 0^") « -► ir-K°)/B(B° -► vr+K-) « 1.2 
is not consistent with the current experimental value of 2.14 ± 0.63 . 

• Since the B ^ b\K decays receive sizable annihilation contributions, their rates are sensitive 
to the interference between penguin and annihilation terms. The measurement of B(B° — > 
bfK~) implies a destructive interference which in turn indicates that the form factors for 
B — > 61 and P — > ai transitions must be of opposite signs. 

• The central values of the branching ratios for the penguin-dominated modes (1270)7r + and 
K-f (1400)7r + predicted by QCD factorization are too small compared to experiment. Just as 
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the case of B — > K*ir decays, sizable power corrections such as weak annihilation are needed 
to account for the observed Kyk rates. The current measurement of B — > ivT-f (1400)-7r + 
seems to favor a K\a — K\b mixing angle of —37° over —58°. 

• The decays B — > i^iK with K\ = ifi(1270) and Ki(1400) are in general quite suppressed, of 
order KT 7 ~ 1(T 8 , except for W -► 7^(1270)^° which can have a branching ratio of order 
2.3 x 10~ 6 . The decay modes K±K + and K~ are of particular interest as they are the 
only AP modes that proceed only through weak annihilation. 

• Time-dependent CP asymmetries in the decays B° — > afn T and bfn T are studied. For the 
former, the mixing-induced parameter S is found to be negative and the effective unitarity 
angle a e g deviates from experiment by around 2a. The discrepancy discrepancy between 
theory and experiment may be resolved by having a larger 7 > 80°. Further precise mea- 
surements are needed to clarify the discrepancy. 

• Branching ratios for the decays B — > fin, fiK, h\ir and h\K with fi = fi (1285), fi (1420) 
and hi = /ii(1170), /ii(1380) are generally of order 10 -6 except for the color-suppressed /i7r 
and hiir modes which are suppressed by one to two orders of magnitude. Measurements 
of the ratios B(B~ -» ^(1380)71") /B(B~ /n(1170)7r _ ) and B(B /1 (1420)^)^(5 
/i(1285)iir) will help determine the mixing angles 0ip 1 and 6z Pl , respectively. 
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APPENDIX A: DECAY AMPLITUDES 



For simplicity, here we do not explicitly show the arguments, M\ and M2, of a? and 0f coeffi- 
cients. The order of the arguments of the oP i {M\M'i) and 0?{M\M-i) is consistent with the order 
of the arguments of the x^ BMl ' M2 \ where 



•f(M 1 M 2 ) = zilgMilMi^. 





Within the framework of QCD factorization [371 ]. the B — > AP decay amplitudes read 

1 



y/2A 



B-^a°Ti- 



V2A 



= % E xi P d) 

V 2 p=M,C 
+ V («1 + 02 

= % E ^ 

V ^ p=u,c 

•^50^0-71+ = X! A p d) 



+ 



^0^+^- = -4 J! A p d) 



^2 p=M,C 



+ 



SpuPi+f% + 



-2ABO^ n o a o = —7= X p d) 
V 2 p=u 



p=u,c 



3 

<5pu («2 - 02) - a 4 + 2 a 3,EW + 2 a 4,EW ~ 03 ~ 03,EW 
+ «I + <EW+/f + /f, E w] 

3 1 

<W (a 2 - 02) - a\ + -af iEW + -a^w - /f - /f, 



(Al) 



x (BTr, ai ) 



+ «4 + «4EW + ^3 + 03,EW 



2 

(S7r,ai) 
1 



EW 



^(Bai,7r) 



<5pu «1 + «4 + «^ EW + /3g + /9 4 2^3,EW o^4,EW 



j^(B7r,ai) 



j^(Sai,7r) I 



A , P _i_ P i nP I /9P flP a?P 

o P u + a 4 + a 4 EYV ■+■ P3 -r P4 — — P3, ew ~~ 2 4> EW 



f,EW 



^-(B7r,ai) 



<V» ("2 - 01) ~ « 4 + T^f.EW + ^°4,EW - ^3 - 2 #4 
1 



1 



+ 2^3,EW 2^,EW 



x (Bir,ai) 



+ 



3 1 

5 P u (a 2 - 0i) - «4 + 2«3,ew + 2 a 4,EW -0l~ 2 0l 
^ 2^ 3 ' EW 2 4 ' EW J I' 



(A2) 



for B — > ai7r, 



•^-B-^a~K° 



a/2 ^ p 

» p=u,c 



5p U /?2 + "4 - -04 EW + 03 + 0l EW 



V * p=u, c I 



jf(Bai,^) 



EW 



+ 



r 1 P 

Opu «2 T 2 a 3,EW 



x (SK, ai ) I 
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A B^a+K- = % E A P } <W Q l + "4 + «4,EW + #3 ~ ^3,EW 



p=u,c 



x (Ba u K) 



V2A 



G F 



E A , s) < 



■«4 2 a4 ' EW ~~ ^ 3 2 3,EW 



x (B ai ,K) 



+ 



5pu OL2 + 2°3,EW 



x (BJ?,oi) I 



and i? — > ai-K", 



s/2A 



G F 



B~^f°n- 



E A f | [W«2 + A) + 2al + a\ + -ag >EW 



- irv p 4- /3 P 4- /? p 1 y(- B7r 'A 9 ) 4- x /0 
o a 4,EW + ^3 + P3,EW ^ + V Z 



a 3 2 a 3,EW 



+ 



W«i + ft) + a\ + a* EW + /f + /3 P EW ] , 



(A3) 



G 

2A B°^f°n« = E A ? } j [W«2 " ft) + 2a§ + + ^ EW 



p=u,c 



+ 



M"«2 - A) + a p - §a p EW - I« P EW + # - i/3 p EW - |^ EW lxW?^l, 



M-./o,- = ^§ E Ai s )|[<5 ptt a 2 + 2a p + Ia p EW ]x(^/ 1 9 ) 

+ ^[«5 pu ft + a§ + al - ^a p EW - i< EW + % + /3 P EW ] , 



+ 



<W«1 + /%) + «4 + «4,EW + ^3 + $S,EW 



^V^F = ^ E A i^ [^2 + 2^ + ^ )EW ]x(^ 
1 V ^ p=u,c L 

+ ^[a p + a p - ia p EW - i« p EW + /f - ^f, EW ]x^\ 



X ( B fl> K ) 
q ) 



+ 



for B -► /^tt and B -► /JVK", 



n< p - -n- p 4 /? p - -tf p 1 X^fi'K') I 
a 4 2 4,EW + P3 2 ^ 3 ' EW J [' 



A B—,K°^- 



^ E A ? M 2 + «4 - .<EW + # + /f ,EW 
V^p=u,c Z 



Gf 



(A4) 



-5 E A ? + A) + < + <ew + /S + Pi 

V ^ p=n,c t 



EW 



X (Bn,Ki) 



+ 



S pu a2 + 2 a 3,EW 



X (BKi,ir) 



A 



p=n,c 
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B 



v " p=u,c 



E4 



X {B-k,K 1 ) 



+ 



(A5) 



for B — > if i7T, and 



A 



i;-_ K - K Q - — 7= E ^ ^pufa + °4 ~ 9«4,EW + ^3 +^3,EW X ( ' BKl ' K \ 
V ^ p=u,c 



A B -^ K o K - = % E A f M 2 +«2-^<EW + /?f + /?3 P EwJ^ {M 

v ^ n=u,C 



i p= n ,c I 



B -^KJK 

A 



w-t;i<" . E A p d) 

V ^ p=u,c 



t| E M 1 + # + ^ewF^ 1 ' + 

V * p=-u, c [ 





2 "4,EW 


KK X J 




2 "i.EW 





p p 

a\ - -at 



LEW + 03 + ^4 2^3,EW 2^' EW 



X (BKx,K) 



A'B*-*K<{K ~ ^ 



fBfKxfK 
G F 



1. 



; 4 2 4,EW 



E A p d) 

'=«,C 

fBfKxfK 



3,EW 2 M4 ' EW 



X (BK,Kx) 



h p - -h p 
°4 2 4,EW 



(A6) 



for B — ► Kiif and P — > K\K, where Ap^ = V^V^, = V p b V* s , and the factorizable amplitudes 
X( BA ^ and X(« P ' A ) are defined in Eq. (|2.48p . The decay amplitudes for B — > bin and b\K are 
obtained from B —* ai7r and aiif respectively by replacing a\ —*b\. Likewise, the expressions for 
B — > /ii7r, /iiK decay amplitudes are obtained by setting (/itt — > /ii7r) and — > h\K). 
The coefficients of the flavor operators a p read 

ai(Mi,M 3 ) = oi(Mi,M 2 ), a 2 {M 1 M 2 ) = a 2 (M 1 M 2 ), 
af(MiM 2 ) = ^(MiM2)-af(MiM2), 
af >EW (M!M 2 ) = ag(MiM 2 ) - a?(MiM 2 ), 

[ a^(MiM 2 ) + r^aliM^)- for MiM 2 = AP, 



o?^M x M 2 ) 



a\^{M x M 2 ) 



a\(M x M 2 ) - r^(MxM 2 ); for MiM 2 = PA, 
a? (Mi M 2 ) + r£ag {M x M 2 ) ; for M x M 2 = AP, 
a? (M!M 2 ) - r^af (MM); for MiM 2 = PA, 



(A7) 



where r B and are defined before in Eq. (|3.12|) . 
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